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Abstract 

We consider general reaction diffusion systems posed on rectangular lattices in two or more spatial dimen- 
sions. We show that travelling wave solutions to such systems that propagate in rational directions are 
nonlinearly stable under small perturbations. We employ recently developed techniques involving point-wise 
Green's functions estimates for functional differential equations of mixed type (MFDEs), allowing our results 
to be applied even in situations where comparison principles are not available. 
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1 Introduction 

A classical line of inquiry in applied mathematics concerns the mechanisms by which energy or 
activity is transported across a spatial domain. In this paper we study the dynamics in higher- 
dimensional spatially discrete models such as the Nagumo lattice differential equation (LDE) 

Uij = Ui+ij + + Uij-1 - 4uij +g{uij), G I?, (1.1) 

with prototypical cubic nonlincarity g{u) = u(u — l)(p — u) for some p G (0, 1). We are interested in 
initial data that are close to a planar travelling wave solution 

Uij{t) = ^iai+ j (72+ ct); $(-00) = 0, $(00) = 1. (1.2) 
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The Nagumo equation is a phenomenological model in which two stable equilibria compete for 
dominance in a spatial domain. In modelling contexts one often thinks of these equilibria as repre- 
senting material phases or chemical or biological species. The LDE (1.1) occurs both as a model in 
its own right as well as a standard semi-discretization of the Nagumo PDE. Planar travelling waves 
(1.2) can be thought of as a mechanism of transport in which the fitter species or more energetically 
favourable phase invades the spatial domain, albeit for very specialized initial data. Establishing 
the stability of these waves, as we do here in this paper, shows that this transport mechanism is 
available for an open set of initial conditions in phase space. 

We note that much of the prior work on stability of travelling waves in (one-dimensional) lattice 
reaction-diffusion equations relies heavily on comparison principles. Indeed, one can often leverage 
the monotonicity inherent in the problem by using Lyapunov functions or sub and super-solutions 
that sandwich the travelling wave; see e.g. [TTl [33] . It is only recently that existence and stability 
results have become available for LDEs that do not admit a comparison principle |3l [20l [21] . In order 
to incorporate such LDEs into our analysis, we avoid the use of comparison principles in our stability 
results. Indeed, in the bulk of the paper we study a general class of models that includes systems of 
LDEs such as the discrete Fitzhugh-Nagumo equation, possibly with nonlinear and non-symmetric 
coupling. For concreteness however, we focus exclusively on the Nagumo equation for the remainder 
of this introduction, but restrict our technical discussion to spectral analysis. 

Existence of Planar Fronts 

To motivate (1.1) and discuss the mathematical difficulties which arise in its study, we recall the 
classical Nagumo PDE 

ut = u^x + Uyy + g{u) (1.3) 

in two space dimensions. This PDE features a thresholding nonlinearity that promotes high frequen- 
cies and has to compete with the diffusion operator, which by contrast attenuates high frequencies. 
The resolution of this competition is the formation of travelling waves, i.e. solutions of the form 

u{t,x,y) = <^{xai + ya2 + ct)\ $(-oo) = 0, $(00) = 1 (1.4) 

which propagate with speed — c in the direction a = (o'i,o'2) with a fixed monotone wave profile 
The existence of such waves can be established via phase-plane analysis, since the wave profile 
necessarily satisfies the planar ODE 

c$' = $" +.g($), (1.5) 

in which we have assumed the normalization (T^ + cr| = 1. 

We remark that in moving from the continuum to the lattice, the first main difference already 
appears when formulating the analogue of the travelling wave ODE (1.5). Indeed, the wave profile 
<& : M M featuring in (1.1)-(1.2) necessarily satisfies the system 

c$'(o = <i>(e + ai) + ct>(e-ai) + <i>(e + <72) + <i>(e-r72)-4<i>(o+,g($(c)); , , 

(1.6) 

$(-00) = 0, $(00) = 1. 

Inspection of this system immediately raises two important considerations that must be addressed. 

The first issue is that due to the shifted arguments, the system (1.6) is a mixed-type functional 
differential equation (MFDE) . The theory for MFDEs is both more subtle and less well known than 
the theory for ODEs. Indeed, the equation (1.6) is not well-posed when regarded as an evolution 
equation on its natural state space C'^([— Cmax, CTmax], II^), with (Tmax = max{|(Ti| , |(T2|}- Nevertheless, 
the theory for MFDEs is sufficiently mature [lH [25l [28l [31] for our purposes in this paper. For 
example, it is known that for each bistable g and for each tr S there is a unique c for which (1.6) 
admits a monotone solution [29 . We emphasize that we do not regard the current paper as making 
any significant contribution towards handling MFDEs. 
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The second issue is that the direction of propagation a features prominently in (1.6), in contrast 
to the ODE (1.5) which is completely isotropic. In fact, this a dependence can be quite rough, 
particularly for choices of g at which the two spatially constant equilibria are nearly but not exactly 
energetically equivalent. For such choices of g the waves move slowly enough to feel the lattice 
discreteness as a kind of friction which further slows down the fronts, sometimes to the point of 
pinning them [71 [HI [H]. Moreover, both numerical and analytical work [101 EH 1131 133 points to 
this additional friction being felt anisotropically, more strongly in directions such as cr — (1,0) 
and a — (1,1) for which there is stronger resonance with the lattice and weakening rapidly as the 
resonance weakens. We regard the main contribution of this work as the extension of techniques in 
nonlinear stability theory to handle the peculiarities that arise in anisotropic systems. 

Stability of Travelling Fronts 

A natural question, and a chief motivator of the present study, regards the extent to which the 
slowing of planar fronts facing lattice directions governs the evolution of initial data which are not 
themselves exact travelling waves. Preliminary numerical experiments show compactly supported 
initial data expanding to occupy a square whose sides are aligned with the lattice directions and 
which propagate outwards at the (slow) speed of the planar fronts. This suggests that the anisotropy 
in the wave speeds of planar fronts governs the dynamics for a large set of initial data. Indeed, one 
may conjecture that fronts facing the lattice directions (1,0) and (0, 1) play a distinguished role in 
governing the dynamics of a large class of initial data. 

In this study we take a small step towards understanding this behaviour by establishing the 
multidimensional stability of planar fronts moving through the lattice 1?' . In particular, we show 
that the wave speeds of the planar fronts govern the dynamics of initial data which are themselves 
close to planar fronts. Furthermore, our results suggest that the waves facing directions in strongest 
resonance with the lattice, i.e. a = (1,0) or (0, 1), are more stable in the following sense. Although 
we obtain the same slow t~^/^ rate for the relaxation of nearly linear interfaces to linear interfaces 
regardless of the direction that the wave is facing, we do obtain the relatively fast t~^/'^ rate for the 
rest of the perturbation to fronts facing a lattice direction as compared to the relatively slow t~^^'^ 
that we are able to obtain for perturbations to fronts facing an oblique direction. 

These observations seem to support the conjecture mentioned above. Of course, much work 
remains to be done and a rigorous examination of initial data that is far from a planar front is 
beyond the scope of this study. Such an analysis may nevertheless be regarded as a natural next 
step, especially in view of the fact that non-planar travelling wave solutions to PDEs have been 
actively investigated of late; see e.g. [171 [M] . 

One-Dimensional Stability In order to illustrate the extra complications that the anisotropy 
of (1.1) imposes on our nonlinear stability analysis, we start by discussing the stability of waves 
in one spatial dimension. The linearization of the one-dimensional version of the PDE (1.3) about 
its travelling wave has a simple eigenvalue at zero coming from the translational symmetry in the 
problem. The remainder of the spectrum is contained away from the imaginary axis in the left half- 
plane. That the only neutral mode arises from translation suggests the Ansatz u{t, x) = ^{x -f- ct -|- 
9{t))+v{x, t). The degree of freedom generated by allowing 9 to vary can be used to demand that the 
rapidly decaying part of the perturbation v is orthogonal to the (adjoint) neutral mode, implying 
that the linear part of the evolution for v generates an exponentially stable semigroup. One can 
then set up and close a bootstrapping argument with v decaying exponentially fast and 9 remaining 
small. 

We remark that the argument sketched above can be transferred from the one-dimensional con- 
tinuum M to the one-dimensional lattice Z. The chief hurdle to be overcome is that unlike in the 
continuum, on the lattice travelling waves are no longer equilibria in a moving frame but rather 
periodic modulo the shift. This observation is not new; techniques to overcome it have been de- 
veloped, for example in [8l 1131 116j . Here we proceed along the lines of the program pioneered by 
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Benzoni-Gavage and coworkers [5] in which Green's functions for hnear LDEs are constructed by 
'fiUing in' the space between lattice points, restoring the relative equilibrium point of view. 

Two-Dimensional Stability In moving from one spatial dimension to two, an additional mathe- 
matical difficulty is that the interface is no longer compact. Instead of occupying a point or interval, 
it now occupies a line or strip. One must now be concerned not only with rigid translations of the 
front but with long wave deformations of the interface which decay slowly at the linear level. The 
slow decay of these long wave deformations manifests as a curve of essential spectrum which touches 
the origin. The presence of essential spectrum is a complicating factor in the nonlinear stability 
analysis of travelling waves because the semigroups associated with the linear part of the evolution 
now exhibit only algebraic decay and thus a bootstrapping argument is no longer guaranteed to 
close if (even bounded) quadratic nonlinear terms are present. Nonlinear stability in the presence of 
essential spectrum which touches the imaginary axis is presently an active area of research; see e.g. 

0- 

Essential spectrum is also present in (1.3), but several authors have managed to either skillfully 
manage this complication [5^ or avoid it altogether [STIES]. Nevertheless, in the current paper we 
need to reconsider this hurdle and develop techniques to deal with the essential spectrum. The chief 
difficulty in our setting is that although the spectral curves are no worse than in the continuum 
case, the spectral projections are complicated by the anisotropy of the lattice. To better see this, we 
now give a brief description of how the proof of stability in [26] can be adapted to the case where 
the front faces a lattice direction a = (1,0) and how this proof must be extended to cope with the 
anisotropy of the lattice when the front faces an oblique direction. 

Lattice Directions In the case when the wave is facing the lattice direction a — (1,0), we use 
the Ansatz 

u,j{t) = $(^ + c< + 9j{t)) + v,j{t). (1.7) 

Here the phase of the wave is allowed to vary in the transverse direction, meaning the direction per- 
pendicular to the motion of the wave and parallel to the interface which separates regions dominated 
by the respective equilibria. After substituting such an Ansatz into the evolution (1.1), one obtains 

= [L,tv],, + $'(0(^?j+i + ^j-i - ^0,) - ^),$'(C + 9,) + N,j{v, 9) (1.8) 

with ^ = i -|- ct, together with the identity 

[Lctv]ij ^ Vi+i^j + Vi^ij + Vij+i + Vij-i - Avij + g'{^{i + ct))vij (1.9) 

and the estimate 

N,,{v,9) = O{{9,+,-0,r + {9,_,-0,)')+O{9,{9,+,+0,_,-20,)) 
+0(0,%) +0(4). 

Our choice of notation emphasizes the fact that the operator L^t is non-autonomous. As a conse- 
quence, it is not immediately clear how to leverage information about the spectrum of Let to obtain 
decay rates on the semiflow it generates. The key is to notice that the operator Let is constant- 
coefficient with respect to j. Performing a Fourier transform in this transverse direction decouples 
modes with different frequencies w and reduces (1.9) to 

[Lu,,ctvioj)]i = w,+i(w) + v,-i(w) + e'"w,(a;) 4- e-'"^^^{uJ) - 4w,;(w) + 5'($(i + ct))v^{Lo). (1.11) 

The work of [31|5] essentially says that the spectrum of £^ —cd+L^,o, regarded as an autonomous 
operator from Ty-'^^°°(M, M) L°°{R,R), governs the evolution semiflow of ^^^(w) = L^,ctv{uj) in a 
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fashion similar to the usual spectral mapping theorem. To obtain information regarding the operator 
it suffices to drop the hats in (1.11) and write 

= -cv\^) + v{i + 1) + v{i - 1) - 2(2 - co^^)v{i) + g\^{i))v{i). (1.12) 

At this point, one notes that is the linearization of the one-dimensional problem about $ which 
has received ample attention from investigators. In particular, it is known that is Fredholm with 
index 0, kernel spanned by and range orthogonal to an adjoint eigenfunction ?/', which can be 
normalized to ($',"0) = 1 [ZHl [29]. Moreover, Lq has no spectrum in the closed right half plane 
other than its simple eigenvalue at zero |131 130| . Furthermore, since depends upon lo only by a 
spectral shift, the spectral projection onto the principal eigenvector is w-independent. This point is 
absolutely crucial. 

Indeed, by demanding 

^Wy(i)?A(i + ct) =0 (1.13) 

i 

uniformly in j, we are essentially taking the equivalent of a spectral projection associated to the 
non-autonomous operator L^t onto a region {Re A < — /?} C C. This implies that L^t generates an 
exponentially stable semiflow on \v \ Vij'ip{i + ct) = 0}. Under this restriction, we will see in §5 
that solutions to (1.8) must satisfy the fixed point problem 

v{t) \ ^ ( e„„(i,0) g,g{t,Q) \ ( v{Q) 
9{t) J y 59^(1,0) gee{t,0) J \ 9(0) ^ 

t( Qvv{tM) Gve{t,to) \ ( MJto;v{ta),e{ta)) ' 



Sevit,to) Geeit.h) ) \ Me (to; v (to), 9 (to)) J °- 

In view of discussion above, we expect to have the exponential bounds Gv*{t,to) — 0{e~^^*'~^°''), 
with * — v,9. In addition, we expect Qe0{t,tQ) to behave similarly to a heat kernel, on account of 
the expansion 

9j = 9j+i + 9j^i - 29 J + 0{v) + 0{9^) (1.15) 
one obtains after imposing (1.13) upon (1.8). The nonlinear terms can be roughljj^ bounded as 

U,{t;v,e)^0{v'') + 0{v9) + 0{{9''f)+0{9e''''), (1.16) 

again with * = w, 9. Here we have 6'* — {9j+i — 9j,6j_i — 9j) and 9j'^ = 9j^i + 9j^i — 29j, while 
multiplication of sequences should be performed pointwise and componentwise. 
In view of the integral inequality 

t 

e"''\l + to)-^dto<C{l + tr^, (1.17) 







we expect that v will decay at the slowest algebraic rate appearing among the terms Afv and Afe- 
Using the heuristic that 9 decays like a one-dimensional heat kernel, which means its fc"' difference 
decays like t we see that the slowest decay in AC comes from the terms (6*^)^ and 06'^*, which 

behave as 

Assuming now that indeed v ^ t^'^/^, one can easily show that 9 ~ t~^^^ by substituting a = 1/4 
and /3 — 3/2 into the integral identity 



il + t~ io)""(l + tor^dto < C(l + t)--nin{"+/9-i.",/9}; when a,/3 ^ 1, (1.18) 
which can be found in [12, Lem. 3.2]. 







^ We are deliberately suppressing all mention of the precise sequence spaces we are referring to. 
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This rough sketch of the a — (1,0) case foUows very closely and matches exactly the results 
obtained there in the case of two dimensions. In fact, in this special case the LDE stability analysis 
is easier than the PDE stability analysis simply because all of the operators are bounded and higher- 
order Sobolev norms are not needed. We remark that this analysis is also reminiscent of the work 
in [5] in which fronts facing the horizontal direction ct = (1, 0, . . . , 0) are shown to be stable in M" 
with n > 4 for the nonlocal equation ut — J * u + g{u) with J appropriately elliptic. 

Oblique Directions The chief difficulty in our work arises when we consider waves that are not 
facing the lattice directions. In order to gain some insight, let us consider a wave facing a general 
direction (cti, 0-2) G with a\ ^ (t| and aia2 7^ 0. We seek an Ansatz which decomposes the general 
perturbation into a slowly decaying perturbation associated with deformations of the interface and 
a rapidly decaying perturbation which accounts for everything else. That is, we write 

Uij{t) ^ ^{n + ct + di{t)) + Vni{t); n = iai+ja2, l^j(Ji-icr2- (1-19) 

Upon substituting this Ansatz into the evolution (1.1) and using the wave profile equation (1.6), we 
obtain 

Vnl = [Lctv]nl + [M,t0]nl + Nni{v, 0) + Oi^' {n + ct + Oi), (1.20) 

where 

-Av^i+g'mn + ct))vni (1.21) 

and 

[M^tO]ni = ^' {n + (Ji + ct){ei+„^ - Oi) + ^' {n + (72 + ct){ei^^, - Oi) 

(1.22) 

+$'(n - ai + ct){ei_,^ - Oi) + $'(n -G2 + ct){9i+,^ - Bi), 

as well as 

Nni{v,e) = ELiO((0i+.,-e/)2 + (0i-.. -^0') 

+ ELi 0(0^(0,+., - ei)) + eLi o(e,(0z-., - ei)) (1.23) 
+o{eiv,,i) + o(vl^). 

The anisotropy of the lattice manifests itself in several places. Firstly, the term MO is linear in 
the first differences of 0, rather than the second differences. In particular, we expect the nonlinear 
estimates 

N,{t-v,e) ^ o{v'') + 0{ve) + o{{o^f) + oiee"), (1.24) 

now with 

= (^i+cr2 ~ (^h&l-ai — Sl,9l^„^ — 9l,9l+cri — &l) ■ (1-25) 

The term 99^ is particularly dangerous, since we can expect it to behave as t^^, which puts us in 
the f3 = 1 case for which (1.18) fails to hold. 

This particular difficulty is often encountered in the study of stability analysis of travelling waves 
in PDE with conservation law structure such as viscous shocks, e.g. [36]. We make use of the simple, 
yet crucial, identity 

ei{9i+i - 90 = l[9l, - 9f - {9i+, - 9i)\ (1.26) 

which can be regarded as a discrete form of uux — (|u^)_^. The square of the first difference decays 
like and hence is safe. The difference of the squares though only decays like t~^. However, we 

can sum by parts, transferring the difference to the heat kernel and in so doing move from the bad 
estimate ^ 

(1 + i - s)-^/\l + s)-Hs - log(l + t){\ + t)-^'^ (1.27) 
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to the good estimate 

t 

{\ + t- s)-3/4(l + s)-i/2^g _ (1 (1 28) 

In particular, we must carefully keep track of the nonlinear terms of order 0{96'^) to make sure that 
they have the correct form. 

Before we can proceed further, we need to decide on a normalization condition to use for the 
phases 9. The linear operators are now given by 

(1.29) 

+g'{<^{n + ct))vn. 

In particular, the dependence on the frequency uj is now non-trivial. It is true that Cq still has a simple 
eigenvalue at zero with the rest of the spectrum in the left half-plane and thus regular perturbation 
yields a curve of eigenvalues A^j and eigenvectors (f)^. However, because the leading eigenvectors 
depend upon oj, the natural orthogonality condition now lives in the Fourier domain. Indeed, it is 
natural to require that for each frequency w, the spectral projection onto 0^ of the Fourier coefficient 
v{uj) must vanish. On the range of such a projection, one would expect the operator L^t to generate 
an exponentially stable semiflow. 

However, we feel that the frequency mixing inherent in such a projection will lead to prohibitively 
complicated expressions for the nonlinearities that will most likely prevent us from verifying the 
conservation-law structure of the 0{69^) terms. This would thwart any attempts to circumvent the 
/3 = 1 problem in (1.18) through the use of (1.26). For this reason, we continue to use the frequency 
independent normalization 

^ w„,V(« + ct) = (1.30) 

n 

that worked so well in the a = (1,0) case. The price that needs to be paid is that the semiflow 
Gve{t,to) now decays algebraically at the rather slow rate of 0((1 + t — to)~'^^'')- The frequency 
dependence of the eigenfunctions (pi^ is directly to blame for this. 

At first glance, this means that the 0{9v) term in (1.24) behaves as t~^, which of course is 
problematic. By carefully inspecting the part of Qyg that decays at the slowest rate, we are able 
to show that the troublesome portion of v is directly proportional to 0^. Together with a detailed 
analysis of the 0{9v) terms in the nonlinearities, this allows the use of (1-26) to prevent the logarithm 
from appearing in (1.27). 

Overview 

In the following sections we flesh out in great detail the sketch presented above. In §2 we state our 
hypotheses and main theorems. Our theorems are of the form "spectral stability implies nonlinear 
stability." Our hypotheses relate largely to the existence of a travelling wave and spectral assumptions 
on the analog of Cq and C^^. Section 3 is devoted to the development of coordinate systems, both 
the (n, I) coordinates for the lattice and the (w, 9) coordinates for u as well as to the derivation of a 
system of evolution equations in these coordinates and the estimation of the nonlinear terms. Section 
4 is devoted to linear analysis, namely pointwise estimates on the Green's functions associated with 
the {v,9) evolution at the linear level. Section 5 is devoted to nonlinear stability, showing that so 
long as all of the terms on the right hand side decay as they ought for t € [0, T], then they actually 
decay as they ought for t e [0,oo). Finally in §6 we show that our theorem is not vacuous. More 
concretely, we verify our abstract hypotheses for the Nagumo equation (1.1). All of the hypotheses 
with analogs in the continuum case [26 ! we can verify analytically. However we need to make an 
additional hypothesis (HSl)^ to control terms which arise due to the anisotropy of the lattice. 
With the exception of the a = (1,0) (resp. a = (1, 1)) case in which these terms are absent (resp. 
benign), we cannot analytically verify this hypothesis and instead use numerics. We conclude in §7 



7 



by discussing the decisions we made that hmit the scope of our work and suggesting potential areas 
for further research. 

Acknowledgments Hoffman acknowledges support from the NSF (DMS-1108788). Hupkcs ac- 
knowledges support from the Netherlands Organization for Scientific Research (NWO). Van Vleck 
acknowledges support from the NSF (DMS-0812800 and DMS-1 115408). 

2 Main Results 

Consider the nearest-neighbour lattice differential equation 

Uij{t) = f{ui+ij{t),Uij+i{t),Ui-ij{t),Uij-i{t),Uij{t)), (2.1) 

for which Uij{t) G M"^ for all G Z^. For convenience, we introduce the symbol tt^ : f°°(Z^;R'^) 
{R'^)^ for any pair G Z^, which acts as 

7T^jU = {ui+ij,Uij+i,Ui-ij,Uij-i,Uij) e {R'^f. (2.2) 

In particular, these operators produce a vector of neighbouring grid points that arc arranged ac- 
cording to the "-|-"-shaped stencil. The system (2.1) can now be written in the condensed form 

Uij{t)=f{^iju{t)). (2.3) 

The conditions we need on the nonlinearity / are summarized in the following assumption. 

(Hf) The nonUnearity / : (M"^)^ ^ M'' is C^-smooth and there exist two points u± e R'^ with 

fiu±,u±,u±,u±,u±) =0. (2.4) 

We emphasize that u- and u+ do not need to be equal. Let us now pick an arbitrary pair (cti , (72) G Z^ 
with gcd((7i,(72) = 1. Our next assumption states that (2.1) admits a wave solution that travels in 
the direction (0-1,(72) through the lattice. 

(H$) There exists a wave speed c ^ and a wave profile $ G C^(]R, W*) that satisfies the limits 

lim $(^) = u± (2.5) 

and yields a solution to (2.1) upon writing 

Uij {t) = $ {iai + ja2 + ct) . (2.6) 

We note that the pair (c, $) solves the functional differential equation of mixed type (MFDE) 

C$'(0 = /(^(^ + ai), $(^ + (72), $(e - ai), $(e - CJ2), $(0)- (2-7) 

In particular, the C^-continuity mentioned in (H<I>) is automatic upon assuming that $ is merely 
continuous. For conciseness, we rewrite (2.7) in the form 

c$'(0 = /([t$](0), (2.8) 

in which the operator r : C(]R,M'') C(M, {R'^'f) acts as 

M(0 = + + ^2).p{^ - ai),p(C - a2),p(e)) e {^"f. (2.9) 
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A standard approach towards establishing the stabihty of the wave (c, <&) under the nonhnear 
dynamics of the LDE (2.1) is to consider the hnear variational problem 

v,j{t) = Df{[T<i>]{iai +ja2 + ct))n+v{t). (2.10) 

As can be seen, the linear operator on the right hand side of this system does not have constant 
coefficients and hence cannot be diagonalized via Fourier transform. It can however be partially 
diagonalized if one takes the Fourier transform in the direction that is perpendicular to the propa- 
gation of the wave, i.e., upon taking iai + ja2 + ct = constant. As explained in §4, one arrives for 
each transverse spatial frequency uj g [— 7r,7r] at an LDE posed on a one dimensional lattice that is 
parallel to the direction of propagation. This LDE is given by 

Vn{t) = I?/([T$](n + ct)) (^e"'-'^v„+,, (t), e-"^i'^i;„+,, it),e''^^^Vn-a^ (t), e'^^'^Vn-a, (i), ^'„(t)) , 

(2.11) 

in which we have n = iai + j(J2- As explained in detail in [201 §2], there is a close relationship 
between the Green's function for the LDE (2.11) and the linear operators 

: iyi'°°(M,C'^) ^ ^""(M,^^), we[-7r,7r], (2.12) 

that act as 

M(e) = -cp'{0+Df{[Tm))[ru,Pm. (2.13) 
in which the operator r„ : C(M,M'^) C(M, iW^f ) is defined by 

[t.pKO = (e^'^^Xe + ai), e-^^^Xe + ^2), e-'^^Xe - 'Ti), e^'^^Xe - a2),p(0) e (K')'. (2.14) 
Written out in the traditional form used in the landmark paper |28) . we have 

5 

M(^) = -cp'(e) + E Au.AOp{(. + rj), (2.15) 



in which the shifts rj and operators Ai^j are given by 



ri 


= 0-1, 




= e-^-i?i/(<i>(C + ai), $ (e + ^2), $(e - ai), $(e - ^2), $(0), 


f2 


= 0-2, 




= e--i-i?2/(*(e + ^1), $(e + ^T2), $(e - <Ji).. $(e - ^2), HO), 


rs 


= -<^i, 


A^AO 


= er^-^^Dsfm^ + di), $(e + ^2), $(e - ai), $(e - ^2), $(0), 




= -0-2, 




= e+-i-i?4/($(e + ^1), $(e + 'Ta), $(e - tn), *(^ - ^^2), *(e)), 




= 0, 


AuAO 


= Drj{H^ + a,), $(e + ^2), $(e - fTi), $(e - fT2), $(0). 



(2.16) 



The formal adjoints of these operators are written as 

Cl:W^^°°{R,C'^)^ L°°{R,C'^), a;eh7r,7r], (2.17) 

and act as 

5 

[Clqm^cq'iO+Y^Al.iOqi^-r,]. (2.18) 
i=i 

Indeed, the designation of formal adjoint is justified by an easy computation which shows that 

(g(0,[/:X(0>rfe= / ([Cg](0,p(C))rf? (2.19) 

-OO J — OO 



9 



holds for all pairs p,q & W^'°° {M.jC'^), where (•, •) denotes the standard inner product on C^. We 
emphasize that for the rest of the paper we will reserve the symbol (•, •) to denote the standard inner 
product on or M'' with £^ (L^) inner products always written out as sums (integrals). 
In view of the fact that rj G Z holds for all the shifts, an easy computation shows that 

e-27Tie{^uj ~ A)e27rif = C^j - 2'Ki£c - X (2.20) 
for all A € C and ^ G Z, in which the exponential shift operator is defined by 

[e.vm^e'-iviO- (2-21) 

In particular, for all uj € [— tt, tt], the spectrum of is invariant under the operation A i— A + 27ric. 
We now proceed to describe the precise assumptions we need to impose on the spectra of these 
operators. 

Since ^{^) approaches u± as ^ — > cx), we can define the characteristic functions 

5 

A±(z) = -cz + ^[ hm A^.iOV^^ (2.22) 

for each uj e [— tt, tt]. Our first spectral assumption states that the characteristic functions associated 
to ~ A cannot have roots on the imaginary axis if Re A > 0. 

(HS)j,j,g For all uj G [— tt, tt] and A G C that have Re A > 0, we have 

det[A^{w) - XI]^0 (2.23) 

for all G M. 

Since (2.20) implies that we can restrict ImA to a compact interval, the a priori bound on solutions 
1/ G M to (2.23) obtained in [24, Lem. 3.1] shows that (2.23) in fact holds whenever ReA > —/3 for 
some /3 > 0. In particular, in view of Thm. A], (IIS)^,^^ implies that the operators ~ \ are 
Fredholm whenever ReA > — /3. 

Notice that tq = r. It is natural to expect that the operator encodes stability properties of 
the wave (c, $) under perturbations that are constant in the direction transverse to propagation. 
Indeed, our next three spectral assumptions on Cq are sufficient to guarantee the nonlinear stability 
of the wave under such perturbations [53]. These assumptions basically state that Cq has isolated 
simple point spectrum at A G {27ricZ} and no other spectrum in the half-plane ReA > — /3. 

(HSl)p There exists P > Q such that £o - A is invcrtible as a map from W^^°°{^, C^) into L°^{R, C^) 
for all A \ 27ricZ that have Re A > — /3. 

(HS2)p The only nontrivial solutions p G 1^1-°° (M, C^) of Cqp = are p = $' and scalar multiples 
thereof. 

(HS3)p The equation Cop = does not admit a solution p G M^i'°°(M, C''). 

Throughout the remainder of this paper, we will use the shorthand (/) :— 

As we have discussed above, the operator Cq is Fredholm. The continuity of the Fredholm index 
together with (IISl)^ imply that the Fredholm index of Cq must be zero. Using Thm. A] together 
with (IIS2)p, one sees that the kernel of Cq must be one-dimensional. Another application of [251 
Thm. A] now in combination with (IIS3)p shows that the kernel of Cq is spanned by a function 
ip G W^"'^'°°(M, M'*) that can be uniquely fixed by the identity 

/oo 
{m,m)d^^l, (2.24) 
-oo 
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since (j) does not lie in the range of Cq. In addition, both (/)(^) and ip{£,) decay exponentially as 
^ — >■ c». Finally, a short computation along the lines of [201 Lem 2.9] yields the useful identity 



/oo 
„ . ,. - oo 



(2.25) 



which holds for every g M. 

Before we proceed, we introduce for any C G M the exponentially weighted spaces 

BCciRX") = {peC{R,C'')\sup^^^e-i\^\\p{0\<^}, 



{pe Ci(M,C'^) I sup5gRe-';i«l[|p(0| + b'(OI] < oo}. 



(2.26) 



The following perturbation result is a direct consequence of the assumptions (HSl)p-(HS3)p. It basi- 
cally states that there is a branch of simple eigenvalues A^^, for Cuj with corresponding eigenvectors 
(pi^ and adjoint eigenvectors ipuj- These eigenvectors decay exponentially at a rate that can be chosen 
to be independent of uj. In addition, with the exception of simple eigenvalues at A;^ + 27ricZ, the 
spectrum of C^i lies to the left of the Hne Re A = — /3. 

Proposition 2.1. Assume that (Hf), (H^), (HS)^^^ and {HS1)^-{HS3)^ all hold. Then there exists 
a constant < <C 1 together with triplets 



(A^,</)^,V^) e C X Tyi'°°(M,C'^) X W^'^i^.C"), 
defined for each uj G (— <^w, 5^), such that the following hold true. 



(2.27) 



(i) There exists /3 > such that for all uj £ {~Si^,Suj). the operator C{uj) — X is invertible as a 
map from W'^^°° (R,^^) into L°°(M,C''), for all X e C that have ReA > -/? and for which 
X-X^i 27ricZ. 



(ii) The only nontrivial solutions p £ W^'°°{\ 
multiples thereof. 

(Hi) The only nontrivial solutions q G W^'°°{I 
multiples thereof. 



3'') of {C^ — Xu:)p = are p = (j)^ and scalar 
Z'^) of {CIj — X^)q — are q — ipui md scalar 



(iv) For all UJ £ {~Sui,Sui) we have 



(2.28) 



In particular, the equation {C^ — Xi^)v — (f)^ does not admit a solution v £ W^'°°{M., C^). 
(v) There exists rj > such that the maps 



UJ I— > 



A^ 



£ BC_^(M,C^) 



(2.29) 



are all well-defined and C'^ -smooth, 
(vi) We have Xq = 0, 4>o — (f) and ipQ — ip. 



Proof. The key observation is that for sufficiently small rj > 0, the operator Cq is also a Fredholm 
operator with index zero when viewed as a map from i?Cl^(M, C*) into i3C_,,(]R, C^). This can be 
seen as in [4j Lem. 4.7] by observing that the operator 



C^:p^ /:o[cosh(7?-)] V(-) 



(2.30) 



11 



is an operator of the form covered by the theory in '28' that in addition is asymptotically hyperbolic. 
In particular, it is a Fredholm operator mapping BC^ (M, C^) into BCo (E, C*) with a one-dimensional 
kernel spanned by the bounded function ^ i-> cosh.{r}^)<j){^) . In addition, the index of Cq viewed as 
a map between these two spaces is zero, which can be concluded from the spectral flow formula in 
[IHl Thm. C]. 

All the properties stated above now follow from regular perturbation. Indeed, the map uj i— > Cu) 
is analytic and the perturbation Ca — C^j is bounded. More explicitly, consider the range 

7^ = Co{BC\{M.,£'^)) c BC_„(M,C'*) (2.31) 

and define J" : 7^ x C x M ^ BC_^(M, C*) by 

F{p,\,uj)^{C^-\){<i>+p). (2.32) 

Observe that J^(0, 0, 0) = and compute the derivative 

i5p,A-^(0, 0, 0) : (p. A) ^ Cop - A0, (2.33) 

which is injective as a consequence of ^ 72. and (1182)^ and surjective because i?C_^(M, C*) — 
7?,© span{(/)}. The desired curves io i— (A;^, <j)u,) are now a direct consequence of the implicit function 
theorem, upon writing (j)^ = (j) + Pu- The remaining curve a; i— )■ t/j^ can be constructed in a similar 
fashion. □ 

The first of our two final assumptions ensures that the curve w — )■ A(^ touches the origin in a 
quadratic tangency that opens up on the left side of the imaginary axis. The second assumption 
guarantees that the only part of the spectrum of the operators with uj G tt, tt] that can 
potentially lie in the half-plane Re A > —77 are the simple isolated eigenvalues A;^ -f l-nicL defined 
for w « 0. 

(HSl)^ Recalling the curves uj i-> (A^^, <^iJ) defined in Proposition 2.1, we have the Melnikov identity 

r 1 (2.34) 

+2/:°^ (1/^(0, \{{dM^^^ - {dM^^^){dM^^^\ (0) ^ 0. 

(HS2)^ For all w £ [— tt, tt] \ {0} and A e C that have Re A > 0, the operator — A is invertible as a 
map from 1^1-°° (M, C'') into i°°(M,C''). 

Proposition 2.2. Consider the setting of Proposition 2.1. We have [dLj\u\u>=Q G i^ o,nd [d^(t>i^]i^^o 
takes values in iM.'^ . In addition, if (HSl)^ also holds, then there exists a constant k > such that 

ReX^<-Kuj^ (2.35) 

for all UJ e {-6^,Si^). 

Proof. To address the imaginary claims we Taylor expand the quantities Ai^, (jjui and ^jj^j in uj about 
w = in the expression (£^ — A^)0^ = to obtain 

O(W^) = Co(f) + Uj(^Co[du,(t)u;]u=Q + [du;C^~d^X^]u,=0(t)^ 

(2.36) 

Taking the inner product against tp and recalling that Cqi/j — leaves 

[dMu=^^ m)A\d^^^]^=o<t>]{0)d^- (2.37) 
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Inspecting the definition of shows that j[9[^£tj]tj=o maps real-valued functions to real-valued 
functions. This together with the fact that and i\) are real-valued implies that i\d^\J\^=i^ is real. 
In order to obtain (2.35), it suffices to integrate the 0(a;^) term in (2.36) against t/j, yielding 

[95A.].=o = /_°l(^(O,[[9S>C.Uo0](O)rf^ 

r T (2.38) 

□ 

Proposition 2.3. Assume that (Hf), (m), {HS)^^^, {HSl )^-{HS3)p and {HS2)^ all hold and recall 
the constant < <5tj <C 1 from Proposition 2.1. Then there exists a constant /3 > such that for all 
A G C that have Re A > —j3 and all uj e [— tt, tt] that have |w| > S^j, the operator Lu: — X is invertible 
as a map from W^'°° {R,C'^) into L°^(R,C'^). 

Proof. This follows from the continuous dependence of Ccu on w, the fact that resolvent sets are 
open and the identity (2.20) which shows we only have to consider spectral parameters A in compact 
sets. □ 

We are now almost ready to state our main results. First, we introduce for any ^? e M the function 

X{'&) e e°^{Z;W^) defined by 

{i + d)). (2.39) 

eel 

This can be seen as the residual of the M'^-valucd function — i [9(^00; ]cj=o after projecting out the part 
along + ■&). 

Theorem 2.4. Consider the IDE (2.1) and suppose that (Hf), (m), (HS)^^^, {HSl)p-{HS3)^ and 

{HS1)^-{HS2)^ are all satisfied. Pick any a < | and 1 < p < oo. Then there exist constants S > 
and C > such that for any u° e £°° {Z'^;R'^') that can be bounde^ by 

there exist -smooth functions 6 : [0,oo) £°°{Z;M.) and w : [0,oo) ^°°(Z^;]R'') that satisfy the 
following properties. 

(i) The function 

Uij (t) = $(ia-i + j(T2 + ct + Oi^^^j^^ (t)) 

+Xi<7i+ja2{ct){0icr2-jai + l{t) - (^)) (2.41) 

'^'Wicri+ja2A(T2~j(Ti {t) 

satisfies (2.1) for all t > with the initial condition 

u,j{0) - $(zai + ja2) + t'?.i+j.2..-2-j^i • (2-42) 

(ii) For all t > and all I Cz 1j we have 

Y,{i'in + ct),wni{t)) ^0. (2.43) 



The £P norm is taken with respect to the first index and the £^ norm with respect to the second index. 
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(in) The function 6 satisfies the estimates 




< C(l + t)-i/4||„o 

< C(l + t)-i/2||„o 



£p(z;£i(Z;R<')) ' 



(2.44) 



(iv) The function ^^iff • [OjOo) — ^°°(Z;IR) defined by 



(2.45) 



satisfies the estimates 



ll^diff(t)|| 
ll^diff(t)|| 



£2(Z;R) 



£°°(Z;R) 



< Cil + t)-^/^\\v° 

< C{l + t)-^\\v°\\ 




£p(Z;£l(Z;R'')) ' 



(2.46) 



(v) The function w satisfies the estimates 




< C(l + t)-"\\v°\\ 

\\ \ \ ( 

< C(l+t)-5/4||^0 




(2.47) 



Our final two results cover situations where the wave travels in directions that are preferred by 
the lattice, in which case the decay rates on v(t) can be improved. Notice that we only recover 
the bounds obtained by Kapitula PB* when analyzing waves that travel horizontally or vertically 
and for which the perpendicular discrete diffusion coefhcients are multiples of the identity matrix. 
In fact, in this case we even get slightly faster decay rates, because we can exploit the embedding 
£^{1;R) C ^2(Z;M) c e°°{Z-R), which does not generalize to the function spaces L«(M;M). 

Theorem 2.5. Consider the LDE (2.1) and suppose that (Hf), (m), {HS)^^^, {HSl)p-{HS3)^ and 
{HS1)^-{HS2)^ are all satisfied. Pick any matrix A G M.'^^'^. Suppose that either ((Ti,i72) = (IjO) 
with 

f{ui+i,j,Uij+i,Ui^ij,Uij-i,Uij) = A[uij+i +5r(uj+ij,-Ui_ij,Uij), (2.48) 

or that (di, (72) = (0, 1) with 

= A[ui+ij +Ui-ij] + g{u,_j+i,Uij_i,Uij), (2.49) 

or that (di, 172) = (1, 1) with 

f{ui+i,j,Uij+i,Ui-i^j,Uij-i,Uij) = A[ui+ij +Ui^j+i + g{uij). (2.50) 

Then we have [9(^0c^]ij=o — 0. In addition, the statements in Theorem 2.4 all hold, after replacing 
(v) by 

(v)' The function v = w satisfies the estimates 



In the special case that a is one of the lattice directions (1,0) or (0,1) and moreover the coupling 
in the transverse direction is linear with weights that are equal both across components in the 
vector u and from one neighbour to the other, we obtain stronger decay. While this case is certainly 
restrictive, it also holds in our motivating example (1.1) and so we state it here as a theorem. 




< C(l + t)-5/4||„0 

< C(l + t)-3/2 ILO 



£P(Z;£1(Z;R'*)) ' 
£p(z:^l(Z;R<')) ■ 



(2.51) 
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Theorem 2.6. Consider the IDE (2.1) and suppose that (Hf), (m), {HS}^^^, and {HSl)p-{HS3)p 
are all satisfied. Pick any constant k > 0. Suppose that either ((Ti,(T2) = (1,0) with 

f{ui+i,j,Ui,j+i,Ui-i,j,Ui,j-i,Uij) = + Uij-i] + g(ui+ij,Ui-ij,Uij), (2.52) 

or that (c7i,0-2) = (0, 1) with 

f{ui+i,j,Uij+i,Ui-ij,Uij-i,Uij) = K[ui+ij +Ui-ij] + g{uij+i,Uij-i,Uij). (2.53) 

Then the triplets (A„, (p^, V'w) are defined for all co G [— tt, tt] and we have (p^ = ^ for all to G [— tt, tt] . 
In addition, the statements in Theorem 2.4 all hold, after replacing (v) by 



(v)" The function v = w satisfies the estimates 



3 Coordinate System 

In this section we develop the coordinate system that we use to analyze the planar nonlinear LDE 
(2.1). We start by introducing the new coordinates 

n = i(Ji+ju2, ^ ^ 

(3-1) 

I = ta2-J(Ti. 

The first of these coordinates represents the direction parallel to the propagation of the wave, while 
the second coordinate represents the direction perpendicular to wave motion. In the sequel we refer 
to n as the wave coordinate and / as the transverse coordinate. 

We emphasize that we always have {n, I) e due to our assumption that (ai, (72) € Z^. However, 
the inverse of the transformation (3.1) is given by 



I 



3 = ['Jl + al]-\na2-la,), ^^'^^ 



which means that the pairs (n, I) in the range of the transformation (3.1) represent only a sublattice 
of 1? . Indeed, rewriting the LDE (2.1) in terms of our new coordinates, we obtain the system 

Unl{t) = fKlU{t)), (3.3) 

in which the linear operator tt^; : {W'-f is given by 

^nl^~ (^^n+(7i,i+o-2 ' ^"+o'2i' — fi ' ci i' — "■2 ' <''2,'+fi ) • (3-4) 

As an example, the new coordinates transform the discrete Nagumo LDE (1.1) into the system 

Unl = [^^U]nl +g{Unl{t)), (3.5) 

in which 

[A^u]„i — Un+ai,l+a2 + Un+a2,l-<Ti + Un-cn,l-a2 + Un-a2,l+cri — '^Unl- (3.6) 

To avoid clutter, we also introduce the operator 

tt'' : £°°(Z^ (R'^f) (3.7) 
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that acts as 

[^><«]„,=7r,>. (3.8) 

The operator tt^ should be regarded as a shorthand for generating a stencil of grid points. Note 
that when (cri,cr2) is equal to (0,1) or (1,0), we have tt^ = 7r+ as introduced in §2. Lifting the 
nonlinearity / in the obvious fashion to act on iW^)^) instead of (M'')^, we can now rephrase 

(3.3) in the compact form 

u{t) = fi-K^'uit)). (3.9) 

In the sequel, we often refer to (3.3) as an equation in local form, while (3.9) is called an equation 
in global form. 

Let us now consider any bounded function C, : M — >■ M*, where we will use both s = d and s = 5d. 
For any shift i? € M, we introduce the sequence G ^°°(Z;M*) that is defined by 

[T^C]n = C(n + r9). (3.10) 

In particular, converts a function defined on M into a sequence by evaluating the function at the 
points Z + 'i?. Furthermore, for any 6» e ^°°(Z;M), we write TqC e e°°{Z'^;W) for the planar sequence 
that is defined by 

[TeCU = C{n + 0i). (3.11) 

By interpreting a shift i9 G M as a constant sequence in £°°(Z;K), we note that we can also write 
T^C € ^°°(Z^;R*) to represent the planar sequence 

[XaCU^Cin + i)), (3.12) 

which is constant in the second coordinate. 

Using this notation, we observe that the travelling wave solution Uni{t) = $(n + ct) to (3.9) can 
be written as 

u{t) = r,t$. (3.13) 

In order to understand the stability of this travelling wave, we set out to seek solutions to (3.9) that 
can be written as 

u{t) = T,t+eit)^ + v{t), (3.14) 

which in local form is equivalent to 

Unlit) = ^{n + ct + 9l{t)) + Vnl{t). (3.15) 

Here the function 

61 : [0,oo) £°°(Z;M) (3.16) 

rcipresents a deformation of the travelling wave profile in the transverse direction, while the remainder 
function 

v:[0,<x,)-^£°°{I?;M.'^) (3.17) 

contains any perturbations that are 'transverse' to this deformation in an appropriate sense. Indeed, 
allowing 6 to vary in (Z; R) buys us Z-many degrees of freedom which we can use to restrict v to 
live in a favourable subspace. 

In order to describe the different norms we employ for the function v, we introduce the sequence 
spaces 

A'p,, = F(Z;£«(Z;M'^)), (3.18) 

parametrized by pairs of exponents 1 < p < oo and 1 < g < oo. A sequence v S Xp^q can be 
interpreted as a planar sequence v € (i°°{7?-,W^) by the identification 

Vnl = {Vn)l- (3.19) 
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In particular, the £p norm is taken along the direction of the wave and the norm is taken in the 
transverse direction. 

Before we can describe the relation between and v, we need to introduce a number of operators. 
First of all, for any a? e M and any 1 < p < cjo, we write 

:e{Z;M.'^) (3.20) 

for the bounded linear operator that acts as 

g^v = ^(V'(n + z?),?;„). (3.21) 

This operator is well-defined as a consequence of the fact that the adjoint eigenfunction decays 
at an exponential rate as ^ — ^ ±oo. This decay also allows us to view as a bounded linear map 

Q^:Xp,q^i''{Z;R) (3.22) 

by writing 

In particular, acts only on the direction of propagation of the wave. Associated to Q^, we can 
introduce the operator 

Pa = [Ti}'p]Q^, (3.24) 
which can be seen either as a bounded linear map 

: F(M,M'^) -i> £P(M,M'^) (3.25) 

or as a bounded linear map 

'■ ^p,q ~^ ^p,q (3.26) 

by writing 

= [TMQ^v]i = <j>{n + ^)Y1 mn' + ^),Vnn)- (3.27) 

n'ez 

For convenience, we also introduce the operators 

Q'^:Xp,q^ei{Z;R) (3.28) 
which arise in the course of computing time derivatives of Qct and act as 

[Q>]; = I] (^''(rJ + t-w). (3.29) 

Note that t/j'iS,) also decays exponentially as ^ — )• ±oo, which can be deduced from the MFDE 
satisfied by tp- Finally, we write 

Pi) ■ ^P,q ~^ ^P,q (3.30) 

for the associated operators = \T^(j)\Q'^ that act as 

[P'Mni = WUQ'Mi = <t>{n + + ^"'')- (3.31) 

n'ez 

Notice that (2.25) implies that 

Q^T^(l> = 1 (3.32) 

for all ?9 e K. In particular, we have the identity P| = P^, both in C{t'P {Z;W^)) and C{Xp^q). 
In fact, the operator P^ can be seen as the spectral projection onto the neutral mode for the 
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linearization of (2.1) about T^^ when posed on the subspace of £°°(Z^;IR'*) that contains only 
sequences that are constant in the transverse direction. Here we are interested in the case when 
arbitrary perturbations are considered, in which case zero is no longer an isolated eigenvalue and the 
spectral projection does not exist in a classical sense. Rather the projection can be thought of 
as a crude approximation of what the spectral projection onto the origin wants to be. See [51137) for 
more examples of extending spectral projections through the essential spectrum. Let us emphasize 
here that is not a projection. 

3.1 Equations for v and 6 

We are now ready to describe the relation between v{t) and 0{t). The goal is to eliminate the neutral 
mode in the best way that we can that does not complicate the eventual Fourier synthesis. To this 
end, we introduce for any 1 < p < oo and q e {2, oo} the subspace 

Xpjt) ^{ve Xp,g I Qatv = e £«(Z;M)} (3.33) 

and demand v{t) G Xjj-g{t). As we shall see below, this condition is sufficient to fix the dynamics of 
both 9 and v as long as these functions are both sufficiently small. In particular, the pair {v,9) € 
X i'^{Z;M.) form coordinates for the solution m in a tubular neighbourhood of the wave profile 

Differentiating the Ansatz (3.15), we find that 

Unlit) = {c + 9i{t))<l>'{n + ct + ei{t)) +v,u{t) 

(3.34) 

= f{[T'i>]{n + ct + 9i{t))) +9i{t)(j){n + ct + 9i{t)) +v„i{t), 

which we abbreviate as 

u{t) - + 0{t)T,t+eit)^ + v{t), (3.35) 

using the convention that the product of a sequence 9 G ^°°(Z;]R) with a planar sequence v S Xp^oo 
is given by [^I'Jni = OiVni- Imposing the condition that the Ansatz (3.15) is a true solution to the 
LDE (3.9), we obtain the following evolution equation for v: 

v{t) = /(7r^T,t+e(t)$ + TT^f (0) - - ^(t)T,t+e(t)0. (3.36) 

It is desirable to write the term f{y) — f{x) as f'{x){y — x) + 0{{y — x)^). To this end we 
analyze the difference Tct{Tr^Tg — Tgr)^ which when added to n^v gives us the relevant {x — y) 
in our Taylor expansion. Note that Tg does not commute with tt^ or r because 9 is in general 
/-dependent. This difference is a vector of differences of the wave profile, the first term of which is 
^(C + (Ti + ^i+CTa) — ^(f + 0"! + 6*0: the other terms being similar except with different indices on and 
signs in front of the aj. To lowest order this term factors as the product $'(^ + f7i)(6'/+o-2 ~ Oi)- This 
is the first time that the anisotropy of the lattice complicates our analysis. Firstly, the coefficient 
on {9i-a2 — di) is $'(^ — CTi) ^ <i>'(^ + CTi) and hence the two slowly decaying first differences 
{(^i+cr-,—Si) and {9i-cr2—9i) do not combine to form a the relatively rapidly decaying second difference 
9i+2+Si-2 — '29i (unless of course cti = 0). Secondly, because the coefficients on the first differences do 
not have identical behaviour in the wave direction, this term cannot be projected away by imposing 
an orthogonality condition on v that is uniform in I. 

We now introduce notation that allows us to decompose the vector of differences TctiTT^Tg—Ter)^ 
and more generally the right hand side of (3.36) into linear terms that generate the evolution 
semigroup, nonlinear terms whose estimation requires care, and nonlinear terms which are estimated 
relatively easily. 

We use the notation 6'* e £°°{Z;M.^) to refer to the sequence 

— (^/+<T2 ^ ^;-<Ti — Oii&i-a2 ^ OiiSi+at — 0) G M^. (3.37) 



18 



In all of our applications the 5th coordinate is zero, but we include the 5th coordinate so that 
dimensions match. A short computation shows that 

[TT^Tet+e^- Wr$]„; = [T{e^)^{- + ct + 9i)\n, (3.38) 

in which we have introduced the symbol 

^iP^ (-^(Tl+ai -^(Tl ) -^(T2+a2 -^<72 5 (7l+a3 (Tl ) (T2+a4 (T2 ? ^) ) (3.39) 

for any a = (ai, . . . , a^) G M^. For any n € Z and i? € M, we now introduce the linear operator 
M{n,'d) : {R'^)^ that acts as 

M{n,^)a = [{aiT^„a2T„„a3T_^„aiT_„„0)(l>{-+^)]^ 

= (ai(j){'& + n + ai), a24>{'& + n + CT2), az4>{'& + n - cti), ai4>{'& + n - 0-2), O) . 

(3.40) 

Upon introducing the nonlinear function : Z x M x ^ (M*^)^ given by 



7W(n,t?;a)= [T{a)'^{- + M{n,d)a, (3.41) 
the fact that is bounded yields the estimate 

M{n,d]a) =0{\a\^), a 0, (3.42) 
which holds uniformly for n e Z and 1? e M. We have hence obtained 

[7r><Tet+e$ - T,t+eT^ni = M{n, ct + 61)6^ + M{n, ct + 0i; Of), (3.43) 
which we write in global form as 

TT^'Ta+e^ - T,t+eT^ = M{ct + 6)6" + M{ct + 6; 6% (3.44) 
Thus (3.36) may be rewritten informally as 

v = L[tt''v + Mr + M]+ Ni-K^'v + MO" + M) (3.45) 
where we have introduced some new notation. The linear operator L is given in local form by 

L{n,-&)v = Df{[T^{n + 'd))v, (3.46) 

while the nonlinear J\f, given by 

N'{n,'&;v) =/([r$](n + i?)+t;)-/([T$](n + '!?)) -£>/([r$](n + i?))t;, (3.47) 
enjoys the estimate 

^f{n,^■,v) =0{\vf), v^Q, (3.48) 

which is uniform for n G Z and ?9 G M and follows from the C^-smoothncss of / together with the 
uniform boundedness of We now begin the process of separating the terms which are easy to 
estimate from those that require more care by writing 

v{t) = Lict + 9m7T>'v{t) + M{ct + 9{t))9O{t)]-0{t)T,t+0(t)^ 
+7^l,c^(w(^),^(^)), 

in which the nonlincarity TZi^ct ■ ^°°(Z2;M'') x e°°{Z;R) -J> £°° {Z'^;R'^) acts as 
[ni,ct{v,9)U = Lin,ct + 9,)M{n,ct + 9i;9f) 

(3.50) 

+N{n, ct + 9i; TT^^v + M{n, ct + 9i)9f + M{n, ct + 9i; 9f)) . 
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We now set out to derive an expression for 6i. Differentiating the condition Qctv{t) = shows 
that 

Qctv{t) + cQ'^vit) = e £°°(Z;M). (3.51) 
Substitution of (3.49) yields the system 

mQctTct+e{t)<t> = QcMct + e{t))[7TXv{t) + M{ct + ^(t))^] + cQ'^v 

+Qctni,ct{vit),0{t)). 

Observe that for each / and at each time t, the quantity QctTct+0i(t)4' is merely a number that can 
be divided away to isolate 6i. In view of the identity QctTct(f> = 1, the uniform bound on = 
implied by H$, Hf, (2.7) and (2.5) and the expansion = 1 — e + O(e^), we can pick a small 
5$ > and introduce functions : {—6g,Sg) — > K that arc fixed by the relation 

[QctTct+e,<P]-' - 1 = -OiQcM + Va{9i). (3.53) 
These functions satisfy the bounds 

V^{6i) = Om^), Oi^O, (3.54) 

uniformly for t? € M and can be extended to map ^°°(Z; {—60,60)) into ^°°(Z;]R). Using these func- 
tions to rewrite (3.52), we obtain 

e{t) = QcMct + 9{t))[Tr''vit) + M{ct + 9{t))e^t)] + cQ'^v{t) 

-9QaTacf>'Q,tL{ct + 9{t))[-K''v{t) + M{ct + 9{t))9^t)] - c9Q,tT,t<f>'Q'^Xt) (3-55) 
+Si,ct{v{t),9{t)), 

in which the nonlinearity 5i,ct : f°°(Z2;]R'^) x £°°(Z; {-60,60)) -)• 1°°{Z;R) is defined by 

Si,ctiv,9) = {l-9Q,iT,t(j)' + Vct{9))Qctni,ct{v,9) +cVcti9)Q'aV 

+Vet{9)QctL{ct + 9) [Tr''v + M{ct + 9)9''] . ^^'^^^ 

Notice that wc arc explicitly excluding terms of order 0{9v) and 0{99^) from the nonlinearity Si^ct 
because these terms require special treatment in our nonlinear stability analysis; see §5. 
Plugging (3.55) back into (3.49), we arrive at 



(3.57) 



v{t) = [I- Pct]L{ct + 9{t))[TT'<v{t) + M{ct + 9{t))9^t)] - cPUv{t) 

+n2,ct{v{t),9{t)), 

in which the nonlinearity 7^2,ct : i°°{Z^;R'^) x £~(Z; {-60,60)) i°°{Z^;R'^) acts as 

Ti2,ct{v,e) = ni,ct{v,9) 

+OTa+0it)<PQctTct4>' QcMct + e{t)) [,,''v{t) + M{ct + 9{t))9\t)] 
+c9Tct+e{t)(t>QctTct<t>'Q'ctv{t) 
-Tct+0(t)(pSi,ct{v{t),9{t)) 
-[Tct+0it) - T,t](ljQctL{ct + 9{t)) [-K^'v + M{ct + e{t))e\t)] 
-c[Ta+0(t) - Ta]ct>Q'ctv{t). 

Notice that wc arc absorbing the 0{99'') and 0{6v) terms into the nonlinearity 1^2,^, in contrast to 
our expression (3.55) for 9. 



(3.58) 
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The next step is to expand the Hnear operators L{'&) and M(i?) around •& = ct. In particular, we 
write 

Let = L{ct), Met = M{ct) (3.59) 

and obtain the system 

v{t) = [I- Pct]Lct[TT'<v{t) + M^te^it)] - cPUv{t) 

+Ti3,ct{v{t),e{t)), 

m = QctLct[n''v{t)+Maeo{t)]+cQ',Mt) (3.60) 

+eH3,ctHt),9'>{t)] 
+S3,ct{vit),e{t)). 

Here we have introduced the hnear operators Hs,ct : e°°{1'^,W^) x £°°{Z;R^) £°°{Z;R) that act as 
H3,ct[v,e''] = Qct[9tfi(i?)]^=ct[7r^t; + Mrtr] 

+QctLct[d^M{'&)U=ae^ (3.61) 
-QcM'QctLct[T^''v + MaO^] - cQaTotcp'Q'aV, 

together with the nonhncaritics TZ^^ct ■ x £°°{Z; {-5e, 6e)) {Z'^;R'^) given by 

+ [I - Pct][L{ct + ei) - LicmTi-^v + M{ct + 9)9^] (3.52) 
+ [/ - Pa]La[M{ct + 9)- M(cf)]r 

and the nonhnearities ^s.ct : f°°(Z2;M'') x £°°(Z; {-5e,5e)) -5- £°°(Z;M) defined by 
Si,ct{v,9) = Si^ct{v,6) 

+Qct [Lict + 9)- L{ct) - 9[d^Li^)]i,^,t] + M{ct + 9)6''] 
+QctLct [Mict + 9)- M{ct) - 9[d^M{'d)]i)=,t] (3.63) 
-9QctTct(^'Qct[L{ct + 9)- L{ct)] [k^'v + M{ct + 0)9^] 
-9QctTa(l>'QctLct[M{ct + 9)- M{ct)]9\ 

Our final step consists of rewriting the linear operator H^ ct in such a way that the nonlinear 
analysis in §5 is simplified. In particular, we reuse the symbol r to write 

[TV\nl = {Vn+cri.l,Vn+„^,uVn-ai,UVn-a^,UVnl) € (3.64) 

for any v e £'^{Z?;W^). Notice that unlike tt^ , the operator r does not act on the transverse direction 
leZ. Indeed, for any v e ^^(Z^; W^) we have the identity 

■k''v = tv + t-v'', (3.65) 

in which 

K]ni = {Vn,l+a2 " '^nU Vn,l-ai " Vnl,Vn,l-cT2 " «nh f^n,i+<7i " Vnh O) , (3.66) 

while the dot product is defined in the obvious manner. A short computation shows that 
QctLctTV = Y.nez{'^i^ + ct),Df{[T<i>]{n + ct))[Tv]ni} 

(3.67) 

= -cQ'ctV- 
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In particular, we see that 

Hi,ct[vM = Qct[d^L{'d)]^=ctT^'<v-QctTct(l)'QctLctT-v^. (3.68) 
Putting everything together, we finally arrive at the system 

v{t) = [I- P,t]L,t[7T''v{t) + M,tr (t)] - cP^Mt) 

+nA,ct{v{t),e(t)), 

e{t) = QctLct['K''v{t)+M,tO''{t)\+cQ'^^v{t) 

+9Hi^a[Tvit), r (t)] + Si,et{v{t), 9{t)). 



(3.69) 



Here we have introduced the linear operators i?4,ct : ^°°(Z^ {W^f) x ^°°(Z;]R5) ^ £°°(Z;]R) that act 



as 



+QctLct [a^M(t9)]^=rt^* - QcM'QctLotMae" ^^'"^^^ 

and are bounded uniformly for t G M. In addition, after picking a smooth cut-off function x"^ : 
[0, oo) [0, 1] that has 

X^([0,1]) = 1, x^X[2,^))=0, (3.71) 



we have introduced the nonlinearities 7^4,ct : ^°°(Z2;M'') x £°°(Z;R) given by 

T^4,ct{v,0) = x'^(||w||^oc(z2.R<i) +25^^ ||6i||^oo(z.K))^3,ct(t^,6') (3.72) 

and the nonlinearities Si^ct ■ £°°(Z2;K'') x £~(Z;M) -> i°°{Z;R) defined by 

S4,ct{v,9) = x'=(||u||^oo(z2.R<i) +25g-^ ||6i||^oo(z.K))53,ct(w,6l) 

+eQct[d^L{'d)]^=ct[T ■ v"] - 0QctTct(p'QctLctT ■ v\ 

We emphasize that the introduction of the cut-offs serves as a notational convenience only. In 
particular, they allow us to formulate the following global bounds on the nonlinearities. 

Proposition 3.1. Suppose that (Hf), (H^), (HS)^^^ and {HSl)^-{HS3),p are satisfied and pick 
1 < p < oo. There exists a constant C > such that for any i G M, G Xp^2 and 9 G £^(Z;M) we 
have the bounds 

(3.74) 
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Proof. Consider two arbitrary sequences a.h E £^(Z;M) and write ab E £^(Z;R) for the product 
sequence [ab]i ~ aibi. Using Cauchy-Schwartz and an elementary estimate it is easy to see that 

lk^ll^i(Z;R) - ll'^ll^2(Z;M) ll^ll^2(z.][{) , l|a6||^2(Z;K) < II"II^2(z.r) ||&||^oo(z.r) • (3.76) 

Another elementary estimate shows that 

lh^^ll^i(Z;R) — ll'^ll^^i'z.^K) ||''||^oo(z.R) , II^^^IL2(Z,R) — II'^II^2(z.r) ||(l||^oc (z-r) ||^'||^oo(z.r) • (3.77) 

The exponential decay of (/)(^) and as ^ — )■ ±00 imply the inclusions Tct4> & £p{Z;R'^) and 

Tct<p' e ^^'(Z;M'*). In addition, the short-range interaction operators tt^ and r are bounded as maps 
from Xp,2 into A'^ 2- Finally, all long-range interactions featuring in TZ^.ct and 54, ct come from the 
linear operators Qd and Q'^. These operators are bounded as maps from Xp^2 into £^{Z;R). In 
particular, the bounds (3.74)-(3.75) follow directly from the estimates (3.76)-(3.77) upon inspection 
of the structure of 7?.4, ct and «S4,ct. □ 

We conclude this section by summarizing our discussion and explicitly stating the relation be- 
tween our modified system (3.69) and the original LDE (2.1). 

Proposition 3.2. Suppose that (Hf), (H^), {HS)^^^ and {HSl)^-{HS3)p are satisfied. Pick any 
to eM. and consider any pair of functions v : ^0,00) — )■ £°°{1?;W^) and 9 : \to^oo) — )■ i!°°(Z;M) for 
which QctoV{to) = and for which the inequalities 

lkWII«°°(Z2;R<i) < 5 

hold for all t > to. Then if the pair {v, 6) satisfies (3.69), we have Qctv{t) = for all t>to and the 
function 

{t) (3.79) 

satisfies the LDE (2.1) for all t>to- 
Proof. Since QciP'd = Q'ct^ can compute 

ii[Qctv{t)] = cQ'^tv{t)+Qdv{t) 

= cQ'dv{t) + QctP'cAt) + Qctni,ct{v{t),e{t)) 

= Qcini,,t{v{t),e{t)) (3.80) 

= Qctn2,ct{v{t),e{t)) 

= 0. 

The last equality can be verified directly from (3.58), in view of the restriction \9i{t)\ < 6g for all 
I e Z, which shows that (3.53) is satisfied for all t>to- The fact that the function u{t) solves (2.1) 
follows directly from the computations in this section. □ 

4 Green's functions 

In this section we focus on the linear part of the system (3.69), which is given by 
v{t) = [I- Pct]Lct[n''v{t) + Mrf^ (t)] - cP^Mt) 

0{t) = QctLa[n''v{t)+MdO'>it)]+cQ'dv{t). ^^'^^ 
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Viewing this system as an ODE posed on £^(Z^;M'') x £^{Z:R), it is elear that sohitions are weU- 
defined in forward and backward time since the hnear operator on the right hand side is bounded 
uniformly for f e K. In particular, we can define a Green's function 

= ( feiio) ) ^ (^^;^^) X ^°°(^^«)) (4-2) 

by writing 

(4.3) 

for the unique solution to (4.1) that has the initial condition v{to) = Vt^ & £°° {Z'^;R'^) and 9{to) = 

Fixing 1 < p < oo, our goal is to restrict attention to initial conditions in the smaller spaces 
X^q{to) X £«(Z;M). In particular, we obtain bounds for the Green's function viewed as a mapping 

Git, to) : X^g, (to) X i^' (Z; M) ^ X^^^^ {t) x t^- (Z; M) (4.4) 

with Qi G {1, 2} and q2 & {2, oo}. 

Our first results show that any solution {v, 9) to (4.1) can in fact be seen as the coordinates of a 
solution w to the system 

jw{t) = L^tTT^'wit) (4.5) 

posed on £°° [I? ; W^) . In particular, there is a close connection between the Green's function Q {t, to) € 
£(£°°(Z2;]R<^)) associated to (4.5) and the Green's function g{t,to) for (4.1). 

Lemma 4.1. Suppose that (Hf), (H^) and {HSl) -{HSS)^ are satisfied. Consider anyvtg G I'^iJ?; 

and Ota € i?°°(Z;E) for which QctoVto = € ^°°(Z;M). Then upon writing w{t) for the solution to 
(4-5) with the initial condition 

w{to) ^ vto + 9toTcto(l), (4.6) 

the pair 

vit) = [I- PctMt), e{t) = Q^w{t) (4.7) 

satisfies (4-1) with 

v{to) = vt„ 0{to)=et,. (4.8) 
Proof. By construction, for any t > to we have 

v{t) = w{t) - e{t)Ta(l>. (4.9) 
Turning first to the equation for 6{t), we compute 

e{t) = cQ'^tw{t) + QcMt) 

= cQ'^^w{t) + QaLaT:-^w{t) (4.10) 
= cQ',tv{t) + ce{t)Q'^tT,t4> + QctLct^'^v{t) + QaLctTi'^e{t)Ta<t>. 

It hence suffices to show that for every 6 £ (Z; M) we have 

cOQ'^M + QctLctTT^'entcP = QaLctM^e^ (4.11) 

To see this, note that the first order terms in the Taylor expansion of (3.43) yield 

TT''eTa(i>-eTaT(i) = Me\ (4.12) 
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In particular, we have 

= L.tMO'^ + ceTa<p\ ^^"-^^^ 

where in the last step we used the fact that £o0 = 0. Furthermore, differentiating QctTct<t> = 1 with 
respect to t, we find 

cQ'^tTctcP + cQcM' = 0, (4.14) 

which establishes (4.11). 

It remains to consider the equation satisfied by v{t). A short computation yields 

= L,tT:'<w{t)-P,tL,t[iT''v{t) + M,te^{t)]-cP'^^v{t)-c9{t)Ta(t>' 
= LaT^''v{t)+L,tTi'^e{t)T,t(t>-PctLct[i^''v{t) + M^t0<'{t)] 
-cPUv{t) - ce{t)T,t(i>' , 

which together with (4.13) completes the proof. □ 

Corollary 4.2. Consider the setting of Lemma 4-1- The Green's function associated to the linear 
system (4-1) can be represented as 



Q{t,to) 



Gvv{t,to) Qve{tM) \^( [I - Pct]G{tM)[I ~ Pcto] [I - Pct]G{t,to)T,j„c^ 
G0v{t,to) g0e{t,to) J { Qctg{t,to)[I - Peto] QctS{t,to)TaJ 



(4.16) 



In order to study the Green's function Q{t,to), we exploit the fact that the LDE (4.5) posed 
on decouples into a set of LDEs posed on Z after taking the Fourier transform in the direction 
perpendicular to the propagation of the wave. In particular, let us introduce for any 1 < p < cxd and 

q e {1,2,00} the frequency space 

-%,q=F(Z;L«([7r,7r],M'')). (4.17) 

For any v E Xp^ n Xp^2 and 9 & (Z; E) D £^(Z; R) we can now introduce the sequence v G Xp^2 and 
the function G L^([-7r, tt], M) defined by 

Vn{uj) = J2 ^nie-'''^, e{uj) = J2 ^/e-''". (4.18) 

In a standard fashion, these maps can be extended to maps v and 6^9 from Xp^2 — ^ '^p,2 and 
f{Z;R) L'^{[-'K,Tr],R) respectively that are bounded and invcrtible via 



v„i 



= ^f e''"vn{uj)du, ^i = ^f e'^^eiu). (4.19) 



The following technical result gives useful bounds for frequency dependent convolution operators. 

Lemma 4.3. Consider any set of maps : £°°(Z;K'') -)> £°°(Z;M'^) defined for w G [-7r,7r] that 
can he represented as 

n,no Vno- (4.20) 

Suppose that there exist constants rj > and ?9 G M for which the bounds 

[KM,no < m(u;)e-''l"-"°+''l , w G [tt, tt] (4.21) 
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hold for some function m(-) e L°"([— tt, tt], K). Then there exists C > that depends only on r) such 
that for any 1 < p < oo and any pair 1 < 92 < 91 < oo, the map 

K : Xp^g^ ^ (4.22) 

defined by 

[Kv]{uj) = K^v{ij) (4.23) 

satisfies the estimate 
in which 



(4.25) 



1 _ 1 1 

S 92 91 ' 

Proof. Writing- Mo, „, = llmll / , , , the generalized Holder's inequality shows 

•' ^ ^1'^= " "C[L1l{[-TT,7r],R),L'32([-TT,7r],R))' ^ H J 

that 

< ||m||^,([_^_^]_R) . (4.26) 

Using Minkowski's inequality, we obtain 

< E„oeze-''l"-"''+''l[/:. \m{uj)r M'^"^ (4.27) 

< E„oeze-''l"-"°+''lM,,,,, ||^;„o(-)IL«([-.,.],M<^) • 

If 1 < p < oo, an application of Young's inequality now yields 



A simpler argument can be used for the remaining case p — oo. □ 
Formally taking the Fourier transform of the LDE (4.5) yields the system of uncoupled LDEs 

dtWi{t,Lo) = Lct^';^lB{t,uj), u;e[-7r,7r], (4.29) 
in which we have introduced the map tt^ G C{(p {Z-.W^) JP [Z; {W^f )) that acts as 

[^^^'In = {e'^-'^v^+^^,e-'^^^Vn+a,,e-'^^'^Vn-.,,e"'^'^Vn-.,,Vn). (4.30) 
We write Qu,{t,to) € £(£p(Z;M'')) for the Green's function associated to the Id LDE 

jw{t) = L,tn:^w{t). (4.31) 

The following preliminary result provides some basic information for the Green's function Gw{t,to) 
that allows us to show that (4.5) is well-posed as an evolution on Xp^2- 

Lemma 4.4. Suppose that (Hf), (m), {HS)^^^, {HSl)p-{HS3)p and {HS1)^-{HS2)^ are all satis- 
fied and pick any e > 0. Then the Green's function Qu;{t,to) can be represented as the convolution 

[GAt,toMn = Yl [OUt,to)]n,noWno (4-32) 

and there exist constants C > and rj > such thai we have the bound 

[Gu.{t,to)]n,no < Ce^(*-*o)e-''l"+^*-"°-^*«l (4.33) 
for every t >to and ui € [— tt, tt]. 
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Proof. Without loss of generality, we assume that c > 0. An important observation due to Benzoni- 
Gavage, Huout and Rousset ^ Thm. 4.2] shows that 



— 1 

[Gu{t,to)]n,no=-^ e^(*-*°)G„(n + ct,no + cio,A)dA, (4.34) 

ZTTZ 



f ■J — tTTC 

for any sufficiently large 7 ^ 1, in which the functions Goj{^^ fo^ ^) solve the MFDE 

(/:„-A)G^(.,eo,A) =5(.-eo) (4.35) 
in the sense of distributions. In view of the symmetry (2.20), we see that 

GUt ^0, A + 27Ttc) - e2^*(««-«)G(e, Co, A), (4.36) 

which implies that the integration contour in (4.34) can be shifted to the left as long as Guj{£,,£,o, ■) 
is analytic. The construction of the functions G^^ in [20l Lem. 2.6] shows that on regions where 

— A is invertible, the map (A,a;) i-^ Gtj(C,5ojA) is continuous and analytic in the first variable. 
In particular, our spectral assumptions allow us to pick 7 = e > for all w S [— tt, tt]. 

The bound (4.33) follows from [201 Eq- (2.35)]. The exponent 77 > can be chosen to be inde- 
pendent of io due to the fact that the characteristic equations 

det[A±(z) - e/] = (4.37) 



have no roots in sufficiently narrow vertical strips |Re2:| < 2r]. This is a consequence of (HS) 
together with the a-priori bound on Imz established in [241 Lem. 3.1]. □ 
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Lemma 4.5. Suppose that (Hf), (m), (HS)^^^, {HSl)p-{HS3)^ and {HS1)^-{HS2)^ are all satis- 
fied. Pick any t^ and Wt^ € '^p,2- Write w{t) for the solution to (4-5) with w{to) = Wt„. Then 
for all t >to we have w{t) G Xp,2, together with 

Wniit) = ^ / e''" [GMo)Uno[^A^)UduJ- (4-38) 



Proof. Let us first consider wtg G Xp^i n Xp^2- The Fourier transform ■^(^(a;) is now well-defined 

for ah oj e [-7r,7r] and Wt^ _ < ||wtj|_:^,^ ^. Upon writing C(*)('^) = Gio{t,to)wt„{^), we see by 

Lemma 4.3 and Lemma 4.4 that C(t) G Xp ^o C Xp 2 for all t > to. Upon writing w{t) G Xp^2 for the 
inverse Fourier transform of C(0, we find by construction that w solves (4.5) with w(to) = Wt„. By 
uniqueness of solutions, we now must have w{t) = w{t) for all t > tg. These conclusions can now be 
extended to wt„ G Xp 2 in a standard fashion. □ 

In light of Proposition 2.1, we can define for any t € M and lu € {—5u:,6uj) the operator 

Qcj,ct{uj)v = ^ Vt^(n + ct)vn, (4.39) 

together with the projection Pj^^ct = Tct4>ujQu),ct- These frequency dependent projections can be used 
to isolate the leading order dynamics of (4.31) associated to the eigenvalues This allows us to 
considerably refine the crude estimates in Lemma 4.4. 

Lemma 4.6. Suppose that (Hf), (m), (HS)^^^, {HSl)p-{HS3)p and {HS1)^-{HS2)^ are all satis- 
fied. Then there exist constants 5^ > Q, C > 0, (3 > Q and rj > Q such that for any uj € [— 7r,7r] with 
\uj\ > we have 

[Gu{t,to)]n,no < Ce-'^(*-*°)e-''l"+'=*-"°-'=*°l, (4.40) 
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while for any lu G (— ^lj) we have 

Quit, to) = e-^"('-*«)rct0c.Qc.,cto + ni^^t, tn), (4.41) 

together with the bound 

|[Hi,^(i,to)]„,„o| < Ce-^(*-*°)e-''l"+^*-"°-^*°l. (4.42) 

Proof. Recall the representation (4.34) and pick a sufficiently small Sx > 0. For convenience, write 

n^{{X,oj)eCx{-S^,S^)\0<\X-XJ<Sx}. (4.43) 

We can proceed as in [201 Lem. 2.7] to write 

GM, ^0, A) = EUC, Co, A) + F„(e, fo, A) (4.44) 

for all pairs (A, a;) € fi, in which iJ^j is explicitly given by 

EU^,^o, A) = -(A - X^r^c^UO^U^oT, (4.45) 

while the remainder (A, uj) i— > Hi^{^, ^q, A) is continuous and analytic in the first variable. In addition, 
we claim that for some C > and 77 > the uniform bound 

|i?<.(C,eo,A)| <Ce-''l«-«ol (4.46) 

holds for all {X,Ld) £ il. To see this, we note that the operator L^, appearing in the proof of [20l 
Lem. 2.7], which in our case can be interpreted as a quasi inverse for £^ — A^^, can be constructed 
in such a way that it depends continuously on uj. This is a direct consequence of the fact that the 
direct sum decomposition Eq. (4.7)] associated to £^ — X^ depends continuously on uj in view 
of Proposition 2.1. 

The meromorphic expansion (4.44) allows us to shift the integration contour in (4.34) to the left 
of the imaginary axis, picking up the appropriate residues. Indeed, the constant /3 > appearing in 
the statement comes from the fact that Propositions 2.1-2.3 imply that C^^ — A is Fredholm with 
index zero on the half-plane Re A > —j3 and is invertible on this half plane except for the simple 
poles at A = Atj which we have now accounted for. 

□ 

From now on, we assume without loss of generality that cti 0, noting that this can always be 
realized by flipping the horizontal and vertical directions if necessary. The key is that we can write 

TctMu,cto = T^tC^Qcto - iar'K^" - l]Tct[du(l)]u=nQcto - ia^^[e'^''^ - l]Tet0[a^Q„,,tJ^=o 
+ Qc.(<,to) 

(4.47) 

for any lo € (— f^w, ^i^), where the C^-smoothness of the branch uj (0^, ^j^) g i?C_r,(M, M'*)^ allows 
us to estimate 

|[Q^(t,io)]„,„ol < C^2g-''l"+-*-"o-cto|^ (4 48) 

Note that we are using e*'^^" — 1 to replace the regular OiiS) term in the Taylor expansion. The 
interpretation of the formal symbol as a difference operator in the transverse spatial direction is 
exploited heavily in the nonlinear argument in §5. 

For any t > 0, we now introduce the function rvit G i°°([— tt, 7r],M) that is given by 

, , f e-^-* UJ G (-(5„,(5„), 
mAuj) = { 4.49) 
1 uji{-&^,8^). 

Exploiting the representation (4.47) allows us to obtain the following estimate. 
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Corollary 4.7. Consider the setting of Lemma, 4-6- There exist constants > 0, C > 0, P > 0, 

K > and ?7 > such that for any uj G [— tt, tt] we may write 

gUt,to) = mt-to{^)Ta(t>Qct, - i<Jl\e"^^ - \\mt-MT,t\du,'^\^=QQct^ 

_i^-i[e-i- _ l]mt_t„(u;)Tet0[5^Q^,eto]c.=o (4-50) 

in which we have the bound 

|[??2,a,(i,io)]n,no| < C[w2e-'^'^'(*-*o) +e-'^(*-*'')]e-''l"+'=*-"''-'=*°l. (4.51) 

Throughout the remainder of this section, we exploit the representation (4.50) to obtain de- 
tailed bounds for the components of the Green's function 0{t, to). We start by providing some basic 
preliminary estimates. 

Lemma 4.8. Fix k> and consider the multiplication operators induced by the functions mk,t{oj) = 
l^^jfe g-Kw^t_ y/jg^g exists C > such that the following estimates hold for all integers k G {0, 1, 2}. 

\\'>T^k,t\\L^{[-Tr,Tr],M)^L-m-TT.,TT],M) — C{l + t) 

< C(l + i)-(W/4), 

\Mmi-.,.m^LHi-.,.m ^ c(i + i)-(w/4), (4-52) 



-(fe/2+1/2) 



Proof. These bounds follow from Lemma 4.3, upon estimating ||TOft,(||j;^s([_^ ^] r) with s = oo for the 
first estimate, s = 2 for the second and third estimates and s = 1 for the last estimate. □ 

For convenience, we write A{t,to) : i^{Z;R) — )• £^(Z;]R) for the Fourier multiplication operator 
that acts in frequency space as 

A{t,to) ■■ 9^ mt-to{(^)e{io), CO € [-7r,7r]. (4.53) 

Using the bounds 

1 



Li([- 



(4.54) 



the estimates in Lemma 4.8 with k = now imply that 

||-^(*)*0)||^2(Z;K)->.^2(Z;R) < C, 

l|-^(*)*0)ll£l(Z;R)^£2(Z.R) < C{1 +t - to)~^^^, 

l|-^(*)*0)||£2(Z;R)^£~(Z;R) < C{1 +t — to)~^^^, 

M(i,io)||£i(Z;R)^£~(Z;R) < C{l + t-to)~^^^. 



(4.55) 



We also introduce the shorthand A'^{t,to) for the operator 9 [A{t,to)9]'^ = A{t,to)6'^, which in 
view of Lemma 4.8 with k = 1 has the bounds 

M*(*7*o)ll<?2(z.R)^<>2(z.R5) < C(l + i-io) 

||^^(i,to)||<?l(Z;R)^<?2(Z.R5) < C(l + 1 - fo)"^^'^, 
\\-^^it,to)\\i2(^2;m)^i°°{Z-M^) ^ C{l+t-to)^^^^, 

M'^(*)*0)||fl(Z;R)_>f=o(Z.R5) < C{1 +t - toy^. 
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Finally, for any nnit vector e G M"'' wo write ^e^(t, to) for the operator 9 \-^\e - A'i^, io)^]*) which in 
view of Lemma 4.8 with fc = 2 has the bounds 

Me^(*i ^o)||£2(z.][j-|^£2(z.][j5) < C(l+t-io) 

Me^(*i ^o)||£i(Z;R)^£2(z.R5) < C(l + i - io)~^^'*, 

l|-4e*(t,to)||£2(z.R)^£oo(z.R5) < C(l+i — to) 

Me^(*i ^o)!|£l(Z;R)^£~(z.R5) < C(l+t-io) 

We now obtain the following estimates on the Green's function to) and the associated operators 
5*(t,fo) that are defined by 

e*(^,^o)(2) = [e(^,^o)(2)]^ (4.58) 

Proposition 4.9. Suwose that (Hf), (m), {HS)^^^, {HS1)^-{HS3)^ and {HS1)^-{HS2)^ are all 
satisfied. Pick any 1 < p < oo and introduce the shorthand 

Xct = - Pct]TctMu.=o- (4.59) 

Then there exists C > such that for any t > to, the Green's function Q{t,to) for (4-1) satifies the 

bounds 

ll^™(i,io)ll;t.,..^Af,.. + ||^-e(t,io)-Xctar'e4-^^(i,to)||,2(z.R)^;t„. ^ C{l + t^t„y\ 
WGvvit, to)\\x,,,^x^,, + \\Sve{t, to) - Xct^l'ei ■ A^it, to) |U(2:;R)^Af,,, ^ Cil + t- ^o)-^/^ 
ll^™(i,to)||;t.,..^Ar,.^ + ||^-e(*'^o)-XctfTr'e4-^^(t,to)||^2(z.R)^;t.^_^ < C{1 + t - to)-^/\ 
\\gvv{t,to)\\x,,,^x,,^ + \\Gve{t,to) - Xctar'e4 • A'>{t,to)\\^,^^.^^^.^^^^ < C{l + t- to)-^/^ 

(4.60) 

in which 64 = (0,0,0, 1,0), together with 

\\Qdv{t, to)|| A'p,2^i?2(z.R) + WQeeit: to) - A{t, io)|li?2(z.R)^£2(z.R) < C(l + t- io)~^^^, 
\\S9v{t,to)\\x^^^^g2(^2.js.^ + \\Geo{t,to) - •^it^to)\\ei{Z;m)~^PiZ;R) ^ C(l + 1 - to)~^^* 
||^/et;(t,to)|lA'p,2^^~(Z;R) + \\See{t,to) - -4(t, to)||^2(z.R)^^oo(z;R) < C{l+t- to) 
\\Qev{t,to)\\ A'p,i^^~(Z;R) + ll^e^(*'*o) - -4(t,io)||<>i(z.R)^^oo(2.R) < C{l+t-to) ^. 

In addition, for any t > to the components of the function G^{t, to) have the bounds 

\\Gvv{tito)\\xp_2^Xl^ + ll%(^'*o)||£2(z.R)^.;t5^ < C(l+i-to)"\ 

\\Gvvi't,to)\\xp^^^xl^, + ll^«e(t'*o)IU(Z;M)^.A'52 < C(l + i-to)"^/^, 

l|e.\(t,to)IU^,,^;,5^ + ||&:,(t,to)||,2(z.R)^;,a^ < C(l + ^-^o)-'^/^ ^^'^^^ 

||ee„(t,to)IU^.,^;,5^ + ||e„%(t,to)||,i(Z;M)^;,5^ < C{1 + t - to)-^/\ 

together with 

ll%(*>*0)IUp_2^.«2(Z;M5) + ||^90(t,to) - -^*(t)*0)||^2(Z;K)^.^2(Z;R5) < C{l+t-to)~^, 

ll^^i)(*'*o)||;tp,n-^2(Z;M5) + ll^9e(*>*o) " -^*(t) *o) ||^i(Z;R)^.^2(Z;R5) < C{1 + t - to)~^^'^ , 
ll^^i)(*'*o)IUp_2^.^~(Z;R5) + ll^ee(*'*o) - -^*(t)*o)||^2(Z;R)^.^<»(Z;lR5) < C(l + t - to)~^^^, 

ll^^i;(*'*o)IUp_j^.^<»(z.R5) + WOeeit^to) - -^*(t)*o)||^i(z;R)^.^~(z;M5) ^ C{l + t- to)'^^"^. 



3/4 (4-61) 
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Proof. In view of the representation of 0{t, to) given in Corollary 4.2, the boimcls above can be clcirivcid 
from the expressions for Qu}{t,to) given in Corollary 4.7. To this end, one can use the estimates in 
Lemma's 4.3 and 4.8, together with the facts that QctTct<P = 1, that [/ — Pct]TctQct = and that 
QctoVto =0. □ 

In the sequel, it is convenient to decompose the function v{t) as 

v{t) = w{t) + XcfJi^ei ■ e\t). (4.64) 

As a direct consequence of Proposition 4.9 we can now obtain the following bounds for the linear 
evolution of w{t). 

Corollary 4.10. Consider the setting of Proposition 4-9. Upon writing 

Qw{t,to) = {Qwv{t,to) Qwe{t,to)) (4.65) 
for the 1x2 operator valued matrix with components 

Qwv{t;U)) ■= Gvvityto) - XctC^l^^i- Q9y{t,to), 
Qwe{t;ta) ■= Gv9it,to) - XctC^l^e^ ■ Qgg{t,to), 



(4.66) 



(4.67) 



there exists C > such that the following estimates hold for all t>to, 

\\Gwv{t,to)\\x^^^x^^2 + \\Ow0{t,to)\\p(^z.,gO^Xp,2 - C'(l + t-io)~S 

ll^«..(i,*0)IU^,,^;,^,, + ||a«,e(t,to)||,i(Z;MH;,^,, < C{1 + t - to)-'^\ 

\\Gn.v{t,to)\\xp^^^x,,^ + \\Gw0{t,to)\\e2^^,^^^x,,^ < C(l + t- ^o)-'/^ 
\\G^vv{t,to)\\xp^^x,^ + \\GMt,to)\\eizm^x,,^ < C{l + t-to)-'/^. 

5 Proof of Main Results 

In this section we set up the nonlinear iteration argument that is used to prove our main results. In 
particular, we fix 1 < p < oo, pick an initial condition 

VoeX^^M, eoei\Z;R) (5.1) 

and consider for f > the fixed point problem 



m ) ' ( Z ) +-''»'^<'''»' ( 9(l,)H4|™(f„),e«(io)l ) * 



(5.2) 



Upon writing 

w{t)=v{t)-Xct(TT'ei-e^t), (5.3) 
we observe by linearity that any solution to (5.2) for some t>0 necessarily satisfies the system 

eo(t,o) \{ vo\,rt( go{t,to) \ { \ 

GUt,0) J\9o J GUt,to) J V eito)H4Tv{to),e<'{to)] J ""'^ 

G^t,to) \ ( n^,cto{v{to),0{to)) \ 

(5.4) 
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We write 

^0 = \\va\\x^_^ + ll^'o||£i(Z;R) • (5-5) 

We also introduce the non-decreasing functions 



together with 



and 



Finally, we write 



M,,2(t) = supo<,„<,(l + tof/^ \Hto)\\x^^, , 

My^^{t) ^ supo<t„<4(l + to)M|t;(to)IUp,^ . 

M,o,2W = supo<t„<t(l + to)'/* ll«^(io)IU^5^ , 

M,o^^{t) = supo<to<ti^+to)'^^\\v%to)\\x^ , 



Me,2it) = supo<i„<,,(l + toY^^ ll^^(to)||,2(z.K) , 

M,,ooW = SUPo<i„<,(l+fo)^/'||^^(to)||,=o(z.M), 

M,.,2(<) = SUPo<i„<,(l + to f/^ ll^^^(to)||^2(z.K5) 

M0'>,oo{t) = SUPo<t^<t(l + to)M|6'*(^o)||^oo(z.K5) 

M„,2(t) = supo<,^<,(l + to)»/*^||u;(io)IU^,, 

M^,^{t) = supo<4„<4(l + fof/"lk(MII;t,„- 



(5.6) 



(5.7) 



(5.8) 



M{t) = M„,2(i)+M^,o„(i)+M^o,2(t) + M„o,^(i) 

+Me,2(0 + Me,ooit) + Meo,2{t) + Meo,oo{t) (5.9) 
+M^,2(t) + M^,^(t). 

Proposition 5.1. S'tippose that (Hf), (m), (HS)^^^, {HSl)^-{HS3)j^ and {HS1)^-(HS2)^ are all 

satisfied and pick any 1 < j> < oo . Then there exists a constant C > such that the following holds 
true. Consider any initial condition {vo,9o) £ X^i{0) x £^(Z;E) and suppose that for some T > 
the functions 

v: [0,T] ^ Xp^2, e :[0,T]^ e'^{Z;R) (5.10) 

satisfy the fixed point problem (5.2) for all < t < T . Then the functions {v, 6) satisfy the nonlinear 
LDE (3.69) for all <t <T and we have the bound 

M{T) < CEo + CM^{T) + CM^{T). (5.11) 

In addition, for all sufficiently small e > there exist functions 

V ■.[Q,T + e]^ Xp,2, ^: [0,T + e] ^•^^(Z;M) (5.12) 

that satisfy the fixed point problem (5.2) for all < t < T + e and have v(t) = v{t) and 0{t) = 9{t) 
for allO<t<T. 

Notice that the exponents 9/8 and 5/4 appearing in Myo^q{t) and M^^q(t) are slighly less than the 

exponents 5/4 and 3/2 that one would expect from the linear evolution of the corresponding variables. 
This is related to the fact that, even if we assume a priori that M{t) is bounded, the nonlinear term 

T^A,cto (i'(*o), ^(^o)) decays at a rate of 0{t(^^) in Xp^i and 0(t(^^^^) in Xp^2- These rates are typically 
too slow to close a nonlinear estimate. The same is true of the term 0{tQ)Hi^ct[Tv{t{^) , 9'^v{ta)\, which 
also decays at a rate of O(to^) in -^^(ZjM). The key is that we will exploit the special structure of 
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the term i74.ct to show that the 0(1^^) decay is not a problem, as is the case when deaUng with 
conservation laws. The slow decay of Ti-i^cto is not a problem for the estimates on My^q, Mg_g and 
Mgo q because the decay of the Green's functions Qyy, Qg^ and Ggy is already faster than that of Gve, 
Gee and Gae^ respectively. 



Lemma 5.2. Consider the setting of Proposition 5.1. For any < Iq < T we have the estimates 
||7^4,cto(«(^o),e(^o))||;,^^ < CM2(to)(l + to)-\ 
\\ni,cto{v{ta),0{to))\\^ < CAP(to)(l + to)-^/\ 



together with 



54xto(«(^0),^(^0))||,l.2.R^ < C[M2(to) + Af3(t„)](l + to)-ll/8^ 



|54.cto("(to),^(to))| 



1(Z;R) 
2(Z;1 



< C[APito) + A'Pitom + to)-'/' 



(5.13) 



(5.14) 



l^(io)i?4.cto[™(io),0]|l,i(z^R) < CM2(^o)(l + ^o)-"/^ 
|e(to)i?4xto[™(io),0]||,2(z.R) < CAP{to){l + to)-'/\ 



(5.15) 



Proof. The estimates for TZcto and Sct^ follow directly from Proposition 3.1, while the bounds for 
i/4,cto follow from (3.76). □ 

We will make use of the following integral inequality, a refined version of (1.18), in order to 
estimate terms on the right hand side of the variation of constants formula. It is often the case that 
given a fixed DuHamel term, there is a tradeoff between how good an estimate one can use on the 
semigroup and how good an estimate one can use on the nonlinear term. In one case we have a 
difference that we can put on either term. In other cases changing the norms will make one term 
decay more rapidly and the other decay more slowly. The point of this lemma is that you don't 
have to choose whether to make the semigroup decay rapidly or whether to make the nonlinear term 
decay rapidly. By using a fast estimate on the semigroup for tg S [0,t/2] and a fast estimate on the 
nonlinear term for to e [t/2,t] you can have it both ways. The observation is far from new, but its 
utility persists. 

Lemma 5.3 (Lemma 3.2 in |12| (restated)). Let a,P,^ and S be positive real numbers, with 7 7^ 1 
and S ^ I. Then there is a constant C > such that 



t/2 



il + t-toy{l + toydto + I il + t-to)-%l + toy>'dto<C{l + t) 

t/2 



\- mm{Q+7-l,(5+;g-l,a,/3} 

(5.16) 



Proof. In view of the fact that 7 7^ 1, we may compute 



Jg^'il + t-toyil + toy^dto 



< 
< 
< 



Itr'^jfii+toydto 



(1 

(1 + 5*)""t^[(1+5*) 



(5.17) 



|l-7l 



^-j^ _j_ min{a,Q:+7 — 1} 



Similarly, using 6^1, we compute 

Jlil + t-tQ)-\l + to)-f'dto < {l + lt)'P Jl^{l + t-to)-^dtQ 



< 
< 



l-<5 



(5.18) 



(l + it)-'°''^{^'^+'5-l>. 



□ 
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Corollary 5.4. Consider the setting of Proposition 5.1. Let S{-,-) denote any of the semigroups 
appearing in Proposition 4-9 and Corollary 4. 10. In addition, let J^{to) denote any of the terms 
Tii,cto{'v(t(i),9{to)), Si^cta[v{to),0{tii)), e[ta)Hi^cta[Tv{to),'d\ or 6'(io)-ff4,cto [™(*o), 0] . Assume that 
the expression S(t,tQ)Af{tQ) is well defined for all < to < t. Let ||A/'(to)||(j denote either the 



£^(Z;M) or Xp q norm of Af{to) as appropriate and let \\S{t,to)\\ 



denote the norm of S{t^to) 



viewed as a map from one o/£'^i(Z;M) or Xp ^-^ into one o/£*^(Z;]R) or Xp^q^, as appropriate. 

Pick constants i^>0, Q!>0, /3>0, 7>0 and S > 0. Write 7'^ = 7— 1 and 7"^ — if 

7 = 1. Similarly, write S"^ = S — 1 if S ^ 1 and S'^ — —v if 5 = 1. 

Then there exists a constant C > such that for any Ci > 0, C2 > and t > 0, the estimates 



\\S{t,to)\\,^ 
together with 

\\S{t,to)h- 
imply the bound 



■q < Clil+t-toY 



< Clil + t-to 



\mto)\\i<C2{l + to)-\ 

mto)\\2<c2{i+to)-^, 



0<ta <t 



0<to<t 



S{t,toW{to)dto 



< cc\C2{i + ty 



- min{Q:+7'^,/9+(5'^,a,/3} 



Proof. Apply Lemma 5.3, replacing 7 with 7 = 1 — in case 7 = 1 and treating S similarly. 



(5.19) 
(5.20) 

(5.21) 

□ 



Our applications of Corollary 5.4 will always have a > 6 and /3 > 7, which is particulary useful 
since 7^^ and S'^ do not appear by themselves in the arguments of the minimum appearing in (5.21). 
In addition, we always have 7 > 1 and /3 > |. Notice that this allows us to show that the part of 
M{t) associated to the terms 

/ 5{t,to) 





e{to)H4[Tw{to),Q\ 



dto + /o 



g{t,to) 

Gw{t, to 




ni,cto{v{to),eito)) 

Si,ctoiv{to),Oito)) 



dto 
(5.22) 

appearing in (5.2) and (5.4) can be bounded as desired. The same holds true for the part of M{t) 
associated to the terms □ 



Jo 



9{to)H4Tv{to),9^to)] 



dto. 



(5.23) 



Indeed, the critical case 7 = 1 occurs in (5.22) only when coupling TZ^ with Q^y, Qov, G^^, Gg,, and 
Gwv, which all have 1 ^ q decay rates a that are faster than those of M^^q, Mg^q, Myo,q, Mgo^q and 
Mw^q respectively. Similarly, the semigroups G^g and GwO appearing in (5.23) have 1 — >■ g decay rates 
a that are faster than those of Myo^q and M^.g respectively. Finally, only Myo^^ and My,^ao feature 
the exponent |. Since /? > f and the 2 — >• 00 norms of G^y^ Gyg, Gwv and Gwe all have S = ^, the 
critical case 6 = 1 occurs only when the decay rate /3 of the ||-||2 norm of the nonlinearity is strictly 
larger than the decay rate of the relevant component of M(t). 

The exact same arguments can also be applied for the components of M(t) associated to 

/ Gyv{t,to) Gye{t,to) - Xctcri^e4 ■ A''{t,to) 
Jo Gev{t,to) Gee{t,to) ~ A{t,to) 

\G^gyit,to) G^ggit,to)-A<^{t,to) 

It hence remains to consider the terms 
/ Xct<JT^ei-A''{t,to) 
Jo A{t,to) 
V A''{t,to) 

These components require a more refined analysis and are discussed in our next result. 





9ito)H4Tv{to),9<'{to)] 



dto. 







d{to)H4[TXcto<^T^e4 ■ 9'>{to),9'>{to)] 



dto. 



(5.24) 



(5.25) 



Note that there is some overlap between the terms appearing in (5.22) and (5.23). 
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Lemma 5.5. Consider the setting of Proposition 5.1. There exists a constant C > such that for 

any unit vector e € and t > 0, the quantities 



satisfy the estimates 



^qit) := /q \\A{t,to)9{ta)0^{to)\\i,^^.^5-j dto, 
^leit) ■■=Jo\\A^{t,to)e-0{to)e-{to) 



dto 



(5.26) 



Ioo{t) 



< CM2(t)(l+t)-l/2^ 

< CM2(t)(l+f)-3/4, 

< CM'^{t){l + t)-K 



(5.27) 



Proof. Restricting ourselves for the moment to Iq{t), we make use of (5.16) with exponents given 
by the following table. 



"a 
7a 

at 
lb 
5 



\\A{t,to)y^e-^{t-to + ir'/* 

\\9-{to)\\%^M{t)^{to + l)-'/' 
\\A'>(t,to)\\ii^i2^{t~to + l)-^/^ 
\\0(to)\\%-^M(t)Hto + l)-'/^ 
\\A{t,to)y^i2 C 

noo\\e\\e^Mitnt-to+i)-'/^ 



\\A{t,t 

\nto)\\i ' 

Wit^tojlUi 

mto)\\% - 

\\A{t,to)h-- 



foc - [t^to + l)-^/^ 

^M(t)2(to + l)-3/2 
_,,oo ^ (t-to + l)-l 

M(i)2(to + l)-l/2 
to. ^{t-to + l)-^'^ 

Mit)Ht-to + l)-^/^ 



The last two rows are straightforward to understand. This corresponds to the time interval [t/2, t] 
on which the semigroup is weak regardless of its rate of decay and the accumulation of the nonlinear 
terms is what is dangerous. Our nonlinear terms ||^?6'^||£2 < ||6'j|£2 ||^?^||^>oo decay so rapidly that wc 
can use even the incredibly weak £^ — > i"^ estimates for A. This is because the exponent S describing 
the decay of A appears on the right hand side of (5.16) only in the form 13 + S — 1 and never by 
itself, thus a large (3 can compensate for a small 5. In our case we have (3 + S —1 = 1/4 for 9 = 2 and 
1/2 for g = oo, enough to establish the desired rates of decay for Ig, as long as the [0, t/2] intervals 
cooperate. 

The interval [0,t/2] is more subtle. As per (1.26) we decompose the nonlinear term 96^ into the 
sum of an (a) term, (0*)^ (where we are interpreting the square of a vector to be the vector of the 
squares of the components) and a (b) term (^^)*. The first of these we estimate in the usual way 
using the £^ £^ norm on the semigroup. For the second of these, we make use of the definition of 
A'^ to transfer the o from 6"^ to A, effectively summing by parts. 

To establish the estimates for I* ^ simply replace A everywhere by A'^ and by ^g*, adding 1/2 
to all values of a and 6. Notice that the critical situations 7=1 and 5=1 are never encountered. □ 

Proof of Proposition 5.1. The fact that solutions to the fixed point problem (5.2) arc solutions to 
the LDE follows from direct differentiation and the linear well-posedness result Lemma 4.5. The 
short term existence claims involving the constant e > can be established in a standard fashion 
by applying the contraction mapping principle. Finally, the bound (5.11) has been verified by the 
discussion above. 

□ 

Proof of Theorem 2.4. In view of the bound (5.11), one sees that whenever Eq is sufficiently small 
the inequality M{T) < 3CEo must hold. In particular, if Ea is sufficiently small it is possible to 
construct a solution to the fixed point problem (5.2) that is defined for all t > by repeatedly 
applying Proposition 5.1. The bounds for 6*^^ stated in item (iii) can be established in the same 
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way as the bounds for 0^{t). The remainder of the decay rates given in the statement of the theorem 
follow directly from the definition of M(t). □ 

Proof of Theorem 2. 5. Upon writing 

= (^/+o-i + ^/-<Ti — '^Sl,0l+a2 + ^(-0-2 ~ 26*;, O), (5.28) 

the estimates for the nonlinearities 7?.4,ct and S4,ct that appear in Proposition 3.1 remain valid upon 
replacing all occurences of \\9\\ by \\9\\ ||^**||. In addition, there exists an operator H^^ct such 
that 6H4^ct[0, = ^-f^5, ct^**- Finally, we have [du(t)w]uj=o = and [9wV'a;]uj=o = 0, which means that 
v{t) = iu{t). In particular, upon writing 

Meoo,2(t) = SUPo<t„<t(l+ to)'/' ||^**(io) 11,2(2.^3), 

Me.<>,oo{t) = supo<,„<,(l + to)^/'||^"*(to)||,oo(z;M3) ^^'^^^ 



and replacing all references to My^q by Mgoo^g, all the nonlinear terms that decay at the rate 0(1^^) 
have disappeared. This allows us to replace the exponent 9/8 appearing in the definition of M^^2 by 
5/4 and the exponent 5/4 appearing in the definition of M^^oo by 3/2. 

With these choices, all the nonlinearities decay at a rate of 0(t(^'^/^) in £^(Z;]R) or Xp,i and 

O(to^/^) in £^(Z;M) or Xp^2- This is sufficiently fast to close the nonlinear argument with direct 
crude estimates. □ 

Proof of Theorem 2.6. The claims concerning the triplets {X^^ , (j)^, , tpu) follow directly upon inspec- 
tion of (2.14), which shows that 

C^p = CoP + 2k{cos{uj) - l)p. (5.30) 

This in turn implies = 2K(cos(ti;) — 1), = and V'lu = V-"! which immediately establishes 
the two conditions (HSl)^ and (HS2)^. Corollary 4.2 now imphes that both Gvv{t,to) and Gv0{t,to) 
decay exponentially instead of algebraically. As above in the proof of Theorem 2.5, the nonlinearities 
decay at a rate of 0{tQ '^^^) in £^(Z;M) or Xp^2- Since these nonlinear decay rates are now the only 
restriction on the decay of w, we can replace both of the exponents 9/8 and 5/4 appearing in the 
definitions of Myj^2 and Myj^^o respectively by 7/4. □ 



6 The Nagumo Equation 

We return now to the scalar equation (1.1) which was discussed in the introduction. The goal of this 
section is twofold: firstly we show that the main theorems are not vacuous by exhibiting an example 
to which they apply; and secondly we present numerical studies which highlight the sensitive nature 
of anisotropy in the problem. 
We again write 

iiijit) = [A+uit)]ij + g{uij{ty, p), (6.1) 
in which the nonlinearity g is given by 

g{u;p) = -^{u^-l){u-p), (6.2) 

with p G (—1,1). The parameter p determines the location of the middle zero of the cubic and can 
be regarded as a measure of the difference in the depth of the wells of the potential G{u) = J g{u)du. 
When p is far from zero, the "local free energy" of the equilibria at ±1 are significantly different 
and hence the more energetically favorable equilibrium invades the spatial domain rapidly. In such 
a parameter regime, the lattice discreteness is hardly felt; phenomenologically (1.1) behaves very 
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Fig. 1: Left: plot depicting Mp{9) = —[5t^uj\uj=o evaluated for (6.1) with detuning parameter p and direction 
of propagation (o-i,(T2) = (cos&,sm&). It is a polar plot, so actually plotted are {cos 6 Mp{6), sin 0Mp{9)). 
The computations were performed with 7 = 10^^. Right: c{p) plot for the fixed direction 9 = 0, which 
corresponds to (cri, (T2) = (1,0). The critical value for p at which pinning sets in is given. The computations 
were performed with 7 = 10~®. 

much like its PDE limit. However, when p is close to zero, the two equilibria are comparably stable, 
the more stable of the two invades the spatial domain slowly, and the lattice discreteness is felt 
keenly. Here we take care to retain the p dependence so that we can contrast the small and order 1 
parameter regimes. Let = P*io') denote the largest value of p for which the wave is pinned (c — O). 

We now verify the hypotheses. As a preliminary step, we state a comparison principle for the 
linearization of (1.1) about its traveling wave, represented in {n,l) coordinates as 

Vnl - {Vn+ai.l~a2 + Vn-cn,l+a2 + "Vyi+a-iJ+ai + Vn-a2,l-<Ti - 4w„i + (/{^{n + ci))w„/} = 0. (6.3) 

In case the equality in (6.3) is replaced with the inequality > (<) we say that w is a super- (sub-) 
solution. We make use of the following standard comparison principle. 

Lemma 6.1. Let K <Z iP' denote the lattice spanned by {01,02) and (ai,— (T2) on which (6.3) is 
naturally posed. Let q",r° £ £°°{A,R) with q^i > r'^i for all {n,l) G A and ^ r°. Let q{t) and 
r{t) denote super- and sub-solutions to (6.3) with initial conditions q^ and r'^ respectively. Then 
Qniit) > Tni{t) for all t > 0. Moreover, if q and r are both true solutions, then the inequality is strict, 
i.e. qni{t) > rni{t) for all t>0. 

We now verify our main hypotheses. 

Lemma 6.2. Consider the model equation (1.1) with nonlinearity g — g{-',p) given by (6.2). Let 

a £l? and let p > p^{a). Then the hypotheses (Hf), (m), {HSl)^, {HS2)^, {HS3)^, (HS)^^^ and 
{HS2)^ all hold. 

Proof. In the context of our abstract framework (2.1) we have f{u) = f{u\ p) with 

/(mi,W2,U3,U4,"5;p) = "1 +U2 +^3 + ^4 + g{u^\p), 

which is clearly and vanishes when all of the Um are equal to u± E { — l,p, 1}. Thus (Hf) is 
satisfied. The notion of p^. along with the existence of a smooth wave profile connecting —1 and 1 for 
arbitrary (even irrational) directions for p > p* has been established in hence (H<I>) is satisfied 
as well. To verify (HS)^,^^ we compute 

^ui'^'^) — A = — -I- 2 cos(i/(Ti -I- wcr2) + 2cos(i/(T2 — t^CTi) — 4-1- g'{u±) — A. (6.4) 
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Since g'{u±) < and since the negative constant —4 dominates the cosine terms, it follows that 
ReA < whenever (6.4) vanishes, establishing (HS)ggj,. The hypotheses (HS2)p and (HSS)^ follow 
from (30| Thm. C] in which the positive function cj) is shown to be the principal eigenvector (in the 
sense of Krein-Rutman) for the positive operator {Cq — X)~^ when A is real and large. Moreover if 
a + il3 is an eigenvalue of Cq, then {a + i/S — X)~^ is an eigenvalue of [£o~''^]~^j which necessarily has 
spectral radius equal to its principal eigenvalue A~^ . Taking A — > oo in the inequality | (a+i/S— A)~^ | < 
A~^ gives the inequality a < 0, thus (HSl)^ is established for eigenvalues with positive real parts. 
To finish verifying (HSl)p it suffices to show, due to the 27ric-periodicity of the spectrum, that with 
the exception of zero, Cq has no eigenvalues on the imaginary axis. 

Suppose that w is an eigenvector of Cq with eigenvalue iv. We will show that v is necessarily 
zero by showing that the eigensolution initially lies below the neutral mode 0, yet touches it at some 
positive time, implying via Lemma 6.1 that the eigensolution is the neutral mode and hence that 
iy = 0. 

First note that the Fredholm Theory of differential-difference operators implies that w is expo- 
nentially localized. Also note that the equations (4.34) and (4.35) establish that (6.3) admits a real 
valued solution of the form wj^j(i) = cos(i^t) Ilew{n + ct) — sin(z^t) Imi«(n -I- ct). Finally, we claim 
that there are positive constants Ki, K2 and K3 with K3 > K2 such that the function 

= - ^20(ri + ct)) + K^^in + ct) (6.5) 

is a super-solution for (6.3) (For a proof of this claim in the a = (1,0) case see [131 Lem. 8.3], 
the proof needs little modification for the general cr case; a related approach for general a can 
be found in [Ml Lem. 7.6]). Since i;+(0) ^ Ki + {K^ - K2)(l){n) > Ki and v^iO) = Rew{n) is 
exponentially localized it follows that v''{t) < v'^{t) for all t > 0, where we have rescaled w so that 
II Re w II 00 < Ki. Taking the limit as A: — >■ cx) along the subsequence tk = yields Re?i;(^) < K3(j){^). 
Let C — mi{K \ Rew(^) < K(f){£^)}. Since C(f){n + ct) is a solution which initially lies above the 
eigensolution v" it does so for all time. Since these solutions touch at t = it follows that they 
are identical, hence that ly — 0. This establishes (HSl)^. 

We now establish (HS2)^. In light of (HS)^^^, the hypothesis (HS2)^ can only fail if there is an 
eigenvalue in the closed right half plane for some w ^ 0. Similar to the argument for (HSl)^, 
we shall assume that for some uj, the operator Cuj has an eigenvalue a + ii^ with a > and 
leverage the comparison principle to show that cu — v = 0. Suppose that Ci^w — {a + iv)w. 
As above, it follows from (4.34) and (4.35) that (6.3) admits a real valued solution wj^;(t) = 
e"* {cosiyt -f ijjI) Yiew{n + ct) — sm{vt + ujI) Imw(n -|- ct)). Comparing to v'^ yields that a < and 
that Rew < ^30. Now restrict attention to a = 0, compare v'^ to Ccf) with C defined as above, and 
use Lemma 6.1 with t = ^ to see that Rew = Ky,(j) and hence that u) = v ~ Q. □ 

The hypothesis (HSl)^ is more difficult to verify, except in the case of special directions of 
propagation. 

Lemma 6.3. Suppose cr ~ (1,0) or a ^ (1, 1) '^i^d, p > p*(cr). Then hypothesis (HSl)^ holds for 
model equation (1.1) 

Proof. We can compute 

[[dlcM=o4>] iO = -<^lm + <^i) + 0(e - ^1)) - ^im + <y2) + 0(e - <y2)) < o, (6.6) 

giving us the desired negative sign for the first term J^{tp{Oj [[d^^i.u]uj=o4']{0)d^ (2.34). This 
computation is valid for all a. 

In the special cases a — (1,0) and a = (1, 1), the map cj i— )• is symmetric with respect to uj, 
which implies that [d^Xi^]uj=o and [9w0lj]lj=o both vanish and hence the second term disappears. □ 

Although our condition cannot be verified for discrete Nagumo with cubic nonlinearity in general 
directions; there do exist other nonlinearities, e.g., nonlinearities similar to the cubic- like nonlinearity 
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Fig. 2: Left: plot depicting Mp(9) — —[5t\iJii^=o , us in Figure 1. The computations were performed with 
7 = 10"*. Right: c{p) plot for the fixed direction 9 — ^, which corresponds to {ai,a2) = (1, 1). The critical 
value for p at which pinning sets in is given. The computations were performed with 7 — 10~^. 



in [T], for which exphcit solutions are known for general rational directions of propagation. To make 
statements for the Nagumo equation with cubic nonlinearity and more general a, we turn to numerics. 
Figures 1 and 2 show the dependence of the Melnikov integral on a and p. For the parameters we 
have investigated we see non- vanishing Melnikov integrals so long as p > p*{cr). However, there 
does not appear to be any obvious structure such as monotonicity that we can exploit to prove 
analytically that the Melnikov integral is non- vanishing. 

The computations were performed using the numerical method developed in |23] . In particular, 
we search for wave solutions of the form 

Uj{t) = Mj-ct), lim $(0 = ±1, (6.7) 

5— >±oo 

by numerically solving the functional differential equatiorj^ 

c$'(0 = 7$"(0 + [$(e±cos6l) + $(^±sin6')-4$(e)]+ 5(^(6; P)- (6-8) 

Notice the extra 7$"(C) term. To make use of robust codes to solve two-point boundary value 
differential equations and to avoid numerical issues that arise in the singular limit c — >■ 0, we follow 
the approach used in [T31[53] and add this extra term in (6.8), with 7 = 10^^, 7 — 10^^ or 7 — 10^'* 
(depending on context; see captions). We refer the reader to [251 for numerical and theoretical 
results concerning the limit 7 — >■ 0. These results strongly suggest that the region of propagation 
failure can be accurately determined by using small but practical values of 7 > 0. In particular, the 
numerical results presented in this section were verified by repeating a subset of the calculations 
with different values of 7 where possible. 

To compute the Melnikov constants, one needs to have detailed information on the adjoint 
eigenfunction V', and the rotated first derivatives Ai := —i[di^Xi^]ui=o and (pi := —i[duj4>bj]uj=o- These 
were also computed numerically. In particular, the actual equations used to determine ip are given 
by 

-7^"(0 - cV''(C) = [V^(C±cos0)+V(?±sin0)-4^/.(O] 

+DgmO)m q^ 



Here the it signs mean both terms are evaluated and added together. 
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with boundary conditions ip{—L) — 0, ip{L) = and ip{0) = 1 for some large L > 0. The actual 
equations used to compute and Ai are given by 

±sin(6')$'(^±cos6') ±cos(0)$'(C=Fsin6l) (6.10) 

A'i(e) = 0, 

again with boundary conditions L) — 0, 4>i[L) — and (/)i(0) — 1. 

7 Discussion 

In this paper we have established the stability of planar fronts travelling in rational directions through 
two-dimensional lattices with nearest-neighbour interactions, under reasonable spectral hypotheses. 
Lattices in three or more spatial dimensions can be treated using the exact same techniques. Based 
on the PDE results [26], we even expect faster rates of decay to occur. This could even considerably 
simplify the nonlinear analysis by eliminating the problematic slowly decaying nonlinearities. 

The restriction to nearest-neighorbour interactions is of course purely artificial and many other 
variants could be considered. For example, in the case of the Nagumo LDE (1.1), it would be inter- 
esting to see what the effects of a cross-shaped Laplacian would be on the decay rate of perturbations 
to waves travelling in the direction (1, 1). 

Our restriction to (cri,U2) £ is of a more serious nature. However, we do believe that the 
spirit of our approach can be extended to irrational directions. Indeed, the difficulties that need to 
be overcome seem to be primarily technical in nature. For example, the Fourier transform in the 
transverse direction would not map the sequence space £^(Z; W^) into L^([— tt, tt], M''), but would be a 
linear map from L^(IR; W^) onto itself. As a first consequence, one can no longer use the embeddings 
t{'L;W^) C £'^{Z;M.'^) C f°°(Z;K''), which implies that any bounds on the nonlinear terms would 
rely on the Sobolev embedding £r^(M;M'') C L°°(M;M''). A second consequence is that one can no 
longer restrict the Fourier frequency uj to the compact interval [— 7r,7r]. This forces a sharpening of 
the condition (HS2)^. However, similar complications were succesfully handled in |26) . emboldening 
us in our view that they will not present a fundamental obstacle in the current setting. 

In establishing our spectral hypotheses for (1.1) in §6, we made extensive use of the fact that 
the Nagumo equation is scalar and enjoys a comparison principle. It is natural to consider what 
can be achieved by using the comparison principles directly and avoiding the intermediate " spectral 
stability implies nonlinear stability" theorem. The work [32] uses comparison principles to establish 
multidimensional stability for the PDE ut = Au -I- g{u) and many aspects of the argument can be 
carried over the lattice. However, after some effort, we have as yet been unable to control new terms 
generated by the anisotropy of the lattice except in the cases a — (1,0) and a — (1,1), in which 
case these terms are absent or benign. We note that [5] , which studies multidimensional stability for 
traveling waves in higher dimensional anisotropic nonlocal equations, also restricts attention to the 
a = (1,0) case and its higher dimensional analogs. The advantage of having an approach based on 
comparison principles is that follow-up questions concerning the basin of attraction of the travelling 
wave and the presence of obstacles are much easier to address. 
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